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Empirical  design  rules  have  been  established,  from  simple  model  experiments, 
to  facilitate  the  proper  design  of  preforms.  Application  of  these  rules  to 
complex  components  requires  an  understanding  of  the  behavior  of  porous  materials 
with  respect  to  densif icat ion,  flow  and  fracture  and  the  effects  of  die  design, 
lubrication  and  temperature. 

An  interactive  computer-aided  design  approach  was  developed  to  facilitate  the 
design  and  evaluation  of  preforms.  The  design  was  accomplished  in  two  phases. 
The  first  phase  described  the  part  to  be  forged  together  with  a  schematic  sub¬ 
division  of  the  part  into  regions.  The  second  phase  involved  the  specification 
of  trial  preform  shapes  for  the  part.  The  program  evaluated  the  input  shape 
for  densif ication  and  non-aggravating  metal  flow.  This  was  accomplished  by 
recognizing  characteristic  deformation  modes  in  the  input  shape  for  which 
available  design  rules  have  been  programmed.  By  interactively  manipulating 
the  preform  shape,  a  proper  preform  may  be  obtained  rapidly.  The  program 
determines  the  size  and  mass  distribution  required  in  the  preform  to  achieve 
ful 1  densif ication. 

The  present  program  was  implemented  for  non-axi symmetr ic  parts  with  a  reasonably 
complex  geometry.  The  procedure  developed  was  applied  to  a  weapon  component 
wherein  the  designed  preform  was  trial  forged  and  metal lurg ical ly  evaluated. 

The  success  of  the  forging  encourages  further  study  to  extend  the  accuracy 
and  capabilities  of  the  present  approach  to  other  part  configurations. 
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Technology  Program  and  was  administered  by  the  US  Army  Industrial 
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1.0  PROBLEM  DEFINITION  AND  SCOPE 


1.1  Introduction 

Powder  forging  is  a  hybrid  process  in  which  preforms 

made  by  conventional  press -and-sinter  powder  metallurgy  (P/M) 

techniques  are  hot  forged  in  closed  impression  dies.  The 

process  is  not  new;  although  production  of  structural  parts 

via  the  powder  forging  route  is  very  recent.  The  earliest 

application  of  this  process,  circa  1910,  was  the  development 

of  tungsten  wires  and  filaments  hot  forged  from  pressure 

fl)  * 

compacted,  pre-smtered  tungsten. 

During  the  last  few  years,  the  automotive  industry 
expanded  the  use  of  conventional  press -and-s inter  methods  to 
parts  with  increasingly  larger  dimensions  and  strength.  The 
requirements  of  size  and  strength  quickly  exceeded  the 
capacity  of  the  conventional  process.  Residual  porosity  was 
the  main  cause  for  the  degradation  of  strength  levels.  To 
meet  the  required  property  standards,  the  feasibility  of 
appending  the  hot  forging  process  to  the  standard  P/M  process 
sequence  of  pressing  and  sintering  was  investigated.  This 
resulted  in  a  significant  reduction  in  the  level  of  residual 
porosity  and,  consequently,  improvement  in  the  strength 
levels.  However,  the  dynamic  properties  of  impact  and 


^Parenthetical  references  placed  superior  to  the  line 
of  text  refer  to  the  bibliography. 


2 


fatigue  strength  fell  short  of  those  in  conventional  cast  or 
wrought  materials. 

Considerable  enhancement  of  dynamic  property  levels  were 
obtained  when  significant  amounts  of  plastic  deformation  and 
lateral  flow  were  permitted  during  hot  forging.  Powder 
preform  forging,  then,  gained  attention  as  a  viable  alternate 
route  for  the  production  of  high  strength  parts.  The  process 
successfully  combined  the  advantages  of  conventional  powder 
metallurgy,  viz,  high  production  rates  superior  material 
utilization  as  a  result  of  elimination  of  much  or  all  finish 
machining  ability  to  form  complex  components  with  excellent 
surface  finish  in  one  forging  blow,  and  the  property 
enhancement  of  true  forging. 

The  requirement  of  lateral  metal  flow  during 
densif ication  of  a  porous  compact  introduces  two  adverse 
conditions  that  make  the  powder  preform  process  somewhat 
complex  and  difficult  to  control.  Firstly,  the  porosity  is 
in  the  form  of  voids  which  provide  ample  sites  for  the 
occurrence  of  fracture.  Secondly,  as  in  the  case  of 
conventional  forging,  excessive  material  flow  leads  to 
pronounced  die  wear.  For  a  successful  process,  suitable 
process  parameter  selection  is  required  so  that  metal  flow  is 
controlled  between  the  lower  limit  to  fully  density  the 
porous  preform  for  acceptable  properties  and  the  upper  limit 
for  fracture  and  excessive  die  wear.  Proper  preform  design 
is  the  single  most  important  parameter,  and  this  has  been 


3 


demonstrated  in  numerous  earlier  studies. '2-4) 

Rational  design  of  powder  preforms  for  complex  shapes 
requires  knowledge  of  densif ication,  flow  and  fracture 
behavior  of  the  material,  and  the  effects  of  die  design, 
lubrication  and  temperature  on  this  behavior.  Preform  design 
guidelines  have  been  established  for  commonly  occurring 
axisymmetric  parts. These  guidelines  were  developed  from 
relatively  simple  model  experiments,  and  successfully  applied 
to  specific  parts.  However,  in  the  present  development  of 
powder  forging,  the  design  of  preforms  to  forge  into  complex 
parts  is  done  by  extensive  trial -and- error  effort.  This 
tr ial-and-er ror  effort  is  usually  the  rate  limiting  step  in 
the  development  of  the  powder  forging  process. 

1.2  Problem  Definition 

The  current  research  program  is  directed  towards  a 
rational  design  of  preforms  for  powder  forging,  without  the 
extensive  laboratory  tr ial-and-er ror  effort.  A  Computer 
Aided  Design  (CAD)  approach  is  followed,  with  the  aim  of 
reducing  the  lead  time  from  the  inception  of  the  process  to 
successful  prototype  forging.  To  demonstrate  the  validity 
and  evaluate  the  merit  of  this  approach,  reasonable  goals 
have  been  set.  Specifically,  the  goals  are  to  (i)  develop  a 
general  approach,  as  far  as  possible,  to  input  the  data 
regarding  the  shape  of  the  finished  part  and  material,  (ii) 
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develop  a  preform  design  scheme  for  the  part,  (iii)  design 
the  preform  for  a  specific  part  using  the  CAD  approach,  (iv) 
forge  the  preform  designed  in  step  (iii),  and  (v)  critique 
the  merit  and  limitations  of  this  approach. 

The  success  of  CAD  of  preforms  for  powder  forging 
depends  upon  the  availability  of  quantitative  information  at 
various  phases  of  the  design.  At  the  present  state  of 
development,  this  information  is  scant.  Available 
information,  in  the  form  of  empirical  guidelines,  is 
predominantly  for  axisymmetric  shapes.  The  specific  part  is  a 
weapon  component  and  is  non-axisymmetr ic.  The  preform  for 
this  part  is  to  be  designed  using  the  proposed  CAD  approach 
and  evaluated.  Applicable  results  of  axisymmetric  parts, 
together  with  limited  information  generated  during  the  course 
of  this  program  are  used  to  formulate  the  preform  design 
scheme.  The  accuracy  of  the  design  process,  then,  is 
restricted  by  the  amount  of  data  that  is  presently  available. 
Further  research  is  needed  to  improve  the  design  strategy.  A 
goal  of  the  current  program  is  also  to  explicitly  define  the 
avenues  of  future  activity. 

The  current  work  is  organized  as  follows: 

(i)  Review  of  the  process  parameters  for  the  powder  forging 
process,  and  identify  the  effects  of  each  parameter  on  the 
overall  process. 

(li)  Review  the  state  of  analytical  techniques  pertaining  to 
deformation  of  sintered  powder  compacts. 
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(iii)  Definition  of  a  logical  basis  for  an  interactive 
computer-aided  design  approach.  This  is  followed  by  the 
implementation  of  the  basis. 

(iv)  Actual  design  of  a  preform  for  the  weapon  component, 
together  with  prototype  forging  and  evaluation. 
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2.0  SINTERED  POWDER  FORGING:  OVERVIEW  OF  PROCESS  PARAMETERS 

2.1  Introduction 

The  pictorial  production  route  for  the  manufacture  of 
powder  forgings  is  illustrated  in  Figure  1.  The  pictorial 
manufacturing  route  for  a  typical  conventional  hot  forging  is 
also  included  for  comparison. In  the  conventional  forging 
route,  the  hot  billet  is  usually  subjected  to  a  succession  of 
blows,  in  a  series  of  dies,  to  develop  the  final  shape. 

Subsequent  machining  operations  follow  to  produce  the 
finished  part,  Figure  2.  In  the  powder  forging  route,  the 
loose  powder  is  compacted  to  a  preform  which  is  sintered  and 
hot  forged,  in  one  forging  blow,  in  one  set  of  dies,  to 
produce  the  finished  part  which  may  require  little  or  no 
subsequent  machining.  This  results  in  a  reduction  of  the 
actual  forging  cost  and  great  improvement  in  press 
utilization.  The  major  reasons  why  this  is  possible  are: 

(1)  The  forging  characteristics  of  a  hot  powder  preform 
differ  significantly  from  those  of  a  wrought  billet. 

(2)  .Metal  can  be  located  in  the  preform  where  it  is 
required,  and  hence  material  redistribution  during  forging  is 
minimized. 

(3)  The  absence  of  flash  reduces  the  forging  load  required. 

The  success  of  a  powder  forging  depends  upon  the  proper 

location  of  metal  in  various  regions  of  the  preform.  This 

i 

! 

! 
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Figure  1.  Pictorial  comparison  of  conventional  forging  and  powder  forging  processes • (5) 
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Figure  2.  Comparison  of  number  of  required  machining  operations  for  conventionally  and 
P/H  forged  parts. (5) 
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influences  the  metal  flow  during  the  forging  stage  and  the 
structural  homogeniety  of  the  finished  part.  The  ductility 
of  porous  preforms  is  severely  limited  due  to  porosity,  hence 
fracture  is  a  severe  problen.  Consequently,  preform  design 
is  critical. 

Two  options  for  hot  working  the  porous  preform  are 
available,  namely, 

(i)  the  preform  has  a  form  that  corresponds  closely  to  the 
finished  part.  In  this  case,  hot  working  is  essentially 
simple  compaction  with  very  little  metal  flow.  This  process 
is  termed  "hot  repressing". 

(ii)  the  preform  has  a  very  simple  shape  but  correct  weight 
and  most  of  the  shape  detail  is  developed  during  forging  by 
metal  flow.  This  process  is  of  interest  to  the  current  inves¬ 
tigation  because  it  imparts  -  better  properties  to  the 
finished  part  than  hot  repressing. 

2.2  Mechanical  Properties 

Production  of  structural  components  from  powder,  on  a 
competitive  basis  vis-a-vis  wrought  products,  requires  that 
the  final  material  be  free  of  residual  porosity.  This, 
however,  is  not  a  sufficient  condition  to  guarantee  optimal 
resistance  to  crack  propagation  in  impact  or  cyclic  loading. 
The  mode  of  densif ication  has  an  important  bearing  on  this 
aspect . ^ 
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In  hot  repressing  of  porous  preforms,  there  is  little  or 
no  metal  flow  in  the  lateral  direction  to  the  applied 
pressure.  Densif ication  ensues  as  a  result  of  simple  pore 
flattening.  In  the  absence  of  relative  motion  between  the 
collapsed  pore  surfaces,  oxides  or  other  contaminant  layers 
on  the  internal  surfaces  of  the  pores  mitigate  against  the 
development  of  a  strong  mechanical  bond  across  the  interface. 
This  results  in  poor  dynamic  properties  for  the  finished 
part.  In  upset  forging,  where  lateral  metal  flow  is 
permitted,  the  presence  of  relative  motion  between  collapsing 
pore  surfaces  leads  to  rupturing  of  the  oxide  layers.  This 
promotes  the  possibility  of  a  sound  metallurgical  bond.(7) 

Static  properties,  viz,  tensile  strength  and  ductility, 
are  comparable  to  wrought  products  in  preforms  that  are  fully 
densified  by  either  repressing  or  upset  forging.  The 
enhancement  of  Izod  values  obtained  by  upset  forging  is 
shown,  in  Figures  3  and  4,  for  varying  degrees  of  lateral 
metal  flow.  Like  toughness,  fatigue  strength  is  also 
increased  by  the  presence  of  lateral  metal  flow. 

The  effect  of  initial  preform  density  on  the  forging 
process  is  evident.  If  material  flow  were  the  only 
consideration,  maximum  toughness  is  to  be  expected  at  some 
intermediate  level  of  initial  preform  density.  The  rationale 
is  that  at  low  initial  preform  densities,  much  of  the 
material  flow  serves  only  to  close  up  porosity  without 
significant  shear  and  interparticle  movement.  Conversely,  at 
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Figure  3.  Impact  resistance  of  601AB  aluminum  alloy,  hot  forged  to 
full  density  as  a  function  of  height  strain  (  and  hence 
lateral  flow).  (3) 
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Figure  k.  Impact  resistance  of  A620  steel  alloy,  hot  forged  to  full 
density  as  a  function  of  height  strain. (3) 
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too  high  a  starting  value  of  preform  density,  insufficient 
deformation  is  imparted  to  the  preform,  for  interpar tide 
displacements,  to  promote  integrity  of  bonding  across 
collapsed  pore  surfaces. 

Superimposed  on  metal  flow,  effect(s)  of  phenomena 
intrinsic  to  the  material  are  also  exhibited.  In  the  case  of 
aluminum  alloy  powders,  preforms  with  lower  initial  densities 
have  greater  pore  surface  area  per  unit  volume  and  hence 
higher  levels  of  the  tenacious  oxide  film.  Consequently, 
preforms  of  92%  initial  density  result  in  higher  levels  of 
toughness  than  those  of  82%  initial  density.  By  comparison, 
it  is  the  lower  initial  preform  density  (75%  cf  85% 
theoretical)  that  results  in  optimum  toughness  in  low  alloy 
steel  powder  forgings.  This  is  attributed  to  the  fact  that 
the  lower  density  preform  have  a  significantly  larger  amount 
of  interconnected  porosity;  during  sintering,  the  oxide 
reducing  gas  can  penetrate  to  a  greater  extent  resulting  in  a 
lower  final  oxide  content. (3)  Toughness  achieves  a 
saturation  value  after  about  50%  height  reduction  for  the 
steel  powder,  whereas,  for  the  aluminum  alloy  powders  no 
saturation  behavior  is  exhibited. 

2.3  Workability 

Sintered  powder  preforms  typically  consist  of  10-30% 
porosity  in  the  form  of  voids  at  the  interstices  between 
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powder  particles.  These  voids  act  as  sites  for  crack 
initiation  and  aid  in  crack  propagation  during  deformation. 
Consequently,  the  ductility  of  preforms  is  severely 
restricted . 

In  order  to  assess  the  extent  of  deformation  that  may  be 
imparted  to  a  material  prior  to  fracture,  the  concept  of 
workability  or  formability  was  introduced.  This  concept  was 
first  put  forth  for  sheet  metal  working, ^  and  then 
extended  to  bulk  forming  processes . 

Workability  is  defined  as  the  degree  of  deformation  that 
can  be  achieved  in  a  particular  metal  working  process  without 
fracture.  In  general,  workability  depends  upon  the  local 
conditions  of  stress,  strain,  strain-rate  and  temperature  in 
combination  with  material  characteristics.  The  stress  and 
strain  rate  in  a  material  undergoing  deformation  is  not 
uniform,  but  varies  from  point  to  point.  This,  in  turn,  is 
determined  by  process  parameters  associated  with  die  design, 
preform  geometry  and  lubrication.  Control  of  these  parameters 
is  used  to  produce  conditions  favorable  for  enhanced 
deformation  prior  to  fracture. 

These  concepts  may  be  qualitatively  expressed  as 
Workability  =  f(material)  x  g(process) 
where,  in  this  relationship,  'f'  is  a  function  of  the  basic 
ductility  of  the  material  and  'g'  is  a  function  of  the  stress 
system  imposed  by  the  process,  in  particular,  the  secondary 
tensile  stresses  generated  during  deformation.  Thus,  a 


workability  analysis  would  include  (i)  a  fracture  criterion, 

' f ' ,  expressing  the  stress  and/or  strain  states  at  fracture 
as  a  function  of  the  strain-rate  and  temperature;  and  (ii)  a 
description  of  the  local  stress,  strain,  strain-rate, 
temperature  histories,  'g',  at  potential  fracture  sites  in 
the  material  during  deformation.  The  workability  analysis 
consists  of  comparing  the  local  conditions  existing  in  the 
material,  for  a  given  set  of  process  variables,  to  the 
limiting  conditions  expressed  by  a  fracture  criterion. 

The  determination  of  the  material  function  applicable  to 
a  broad  range  of  stress  and  strain  states  is  critical  to  a 
workability  analysis.  Lee  and  Kuhn^^^  successfully 
determined  such  a  criterion  from  upset  tests  on  circular 
cylinders.  The  upset  test  is  particularly  flexible  as  by 
varying  (i)  the  cylinder  aspect  ratio  (height  to  diameter 
ratio)  and  (ii)  the  lubrication  at  the  die  contact  surfaces, 
the  barrelling  severity  and  hence  the  secondary  tensile 
stresses  developed  on  the  bulge  surface  can  be  varied  over  a 
wide  range.  The  results  are  shown  in  Figure  5.  Strain  paths 
are  obtained  from  grid  measurements  taken  at  the  bulge 
surface  on  an  equitorial  diametral  plane.  The  strain  paths 
are  process  dependent  and  their  terminal  points,  determined 
by  visual  observation  of  fracture,  appear  to  lie  on  a 
straight  line  with  a  slope  of  one-half.  This  linear 
relationship  has  been  demonstrated  for  a  variety  of 
materials,  including  sintered  powder  materials,  both  during 
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STRAINS  AT  FRACTURE 


Experimental  fracture  strains  in  bulk  forging  of  wrought 
1045  steel.  (11) 


Figure  5- 
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hot  and  cold  forging. The  intercept  of  the  line  with  the 
tensile  strain  axis  changes  in  accordance  with  the  inherent 
ductility  of  the  material. 

The  -workability  limit  diagram  , Figure  5,  applies  for 
free  surface  fracture.  The  presence  of  hydrostatic  pressure 
tends  to  retard  the  initiation  of  Eracturc.  It  has  been 
shown^  that  the  workability  diagram  is  raised,  towards 
larger  deformation  to  fracture,  in  the  presence  of 
hydrostatic  pressure. 

2.4  Fracture 

A  variety  of  fractures  and  defects  occur  in  a  sintered 
powder  forging  process.  Their  location  and  occurrence  are,  in 
general,  controlled  by  the  process  parameters.  Typical 
fractures  that  occur,  are  illustrated  in  Figure  6  and,  may  be 
classified  as  follows: 

2.4.1  Free  Surface  Fracture 

This  type  of  fracture  occurs  on  the  expanding  free 
surface  of  a  deforming  preform,  before  contact  with  the  die 
walls.  The  cracks  developed  on  the  bulge  surface  of  an  upset 
cylinder  is  typical  of  this  type  of  fracture.  The  occurrence 
is  attributed  to  the  secondary  tensile  stresses  developed  due 
to  friction  at  the  die/workpiece  interface. 
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2.4.2  Die  Contact  Surface  Frature 

This  class  of  fracture  occurs  at  the  die/workpiece 
interface.  A  common  location  for  the  initiation  is  in  the 
vicinity  of  a  punch  or  die  corner  where  metal  flow  changes 
direction.  They  do  not  usually  propagate  throughout  the 
workpiece,  but  does  severely  detract  from  the  surface  finish. 

Fundamental  studies  on  hub  extrusion  by  Suh(4>  show 
that  a  concentration  of  shear  deformation  combined  with 
tension  or  low  values  of  hydrostatic  pressure  in  the  vicinity 
of  a  corner  lead  to  the  initiation  of  this  fracture. 
Non-uniform  lubrication,  also,  plays  a  role.  The  occurrence 
may  be  suppressed  by  providing  generous  radii  and  uniform 
lubrication  at  these  sites,  as  both  contribute  to  reducing 
the  magnitude  of  tensile  stresses. 


2.4.3  Internal  Fracture 


Unlike  free  surface  and  die  contact  fractures,  internal 
fractures  occur  inside  the  workpiece  and  rarely  penetrate  to 
the  surface.  As  a  result,  they  may  go  undetected  and 
possibly  lead  to  catastrophic  failures  during  service. 

Central  bursts  during  extrusion  are  an  example  of  this  class 
of  fractures. 

The  center  of  double  hub  extrusions  is  a  potential  site 
for  an  internal  fracture.  Experimental  studies  of  Suh^4) 
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indicate  that  the  occurence  of  internal  fracture  correlate 
with  the  geometry  of  the  deformation  zone.  Die  lubrication 
has  little  effect. 

2.4.4  Intrusion  Defects 

Intrusion  defects,  although  not  fractures,  also  occur 
during  forging.  They  manifest  as  a  surface  disturbance,  and 
detract  from  a  uniform  surface  finish.  The  occurrence  of 
these  defects  are  generally  geometry  related,  with 
lubrication  having  little  effect. 

2.5  Fracture  Prevention 

A  comprehensive  study  of  fractures  and  their  prevention 
for  simple  axisymmetric  powder  forgings  is  presented 
elsewhere. ^  As  mentioned  in  the  earlier  sections, 
deformation  zone  geometry,  and  consequently  the  preform 
geometry,  have  an  important  bearing  on  the  likelihood  of 
fracture.  Due  to  the  limited  ductility  of  a  porous  preform, 
control  of  metal  flow  becomes  critical.  The  metal  flow  that 
occurs  under  actual  conditions  is  determined  by  the  preform 
shape  and  its  relation  to  the  shape  of  the  dies.  A  subtle 
change  in  the  preform  shape  can  significantly  alter  the  metal 
flow  and  the  likelihood  of  fracture.  This  fact  is 
illustrated  in  Figure  7. 
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Figure  7. 


Preform  options  for  forging  an  axisymmetric  component 
consisting  of  a  flange  and  hub  with  central  hole. (3) 
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The  component  forged  is  axisymmetric,  consisting  of  a 
hub,  a  flange  and  a  central  hole.  Three  preform  options  may 
be  pursued.  Option  1  is  to  have  the  preform  fill  the  flange, 
so  that  the  hub  is  formed  by  back  extrusion.  Option  2  is  to 
have  the  preform  fill  the  hub,  so  that  the  flange  is  formed 
by  radially  outward  metal  flow.  Finally,  option  3  is  an 
intermediate  shape,  where  metal  flow  is  directed  radially 
inward  toward  the  mandrel  and  up  into  the  hub. 

The  likelihood  of  free  surface  fracture  in  each  of  these 
preform  options  may  be  examined  by  studying  the  strain  states 
at  the  free  surface  along  with  the  workability  diagram. 

Figure  5.  Typical  strain  states  which  occur  are  shown  in 
Figure  7.  In  preform  option  1,  the  strain  path  is  located  in 
the  first  quadrant,  the  biaxial  strain  state  being  developed 
due  to  frictional  constraint  at  the  mandrel  and  the  hub 
surfaces,  and  fracture  is  a  strong  possibility.  In  preform 
option  2,  the  strain  path  is  in  the  second  quadrant,  hence  a 
larger  deformation  to  fracture  may  be  expected.  In  option  3, 
the  strain  state  developed  is  biaxial  compression  and 
fracture  is  remote. 

Thus  from  qualitative  arguments  based  on  observation  of 
simple  model  experiments,  a  workable  preform  geometry  may  be 
inferred.  Experiments  have  verified  that  preform  option  3, 
in  the  above  example,  results  in  a  defect  free  forging, 
whet  ;as  the  other  two  preform  options  exhibit  the  fractures 
as  expected . ^ 
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2.6  Lubrication  and  Temperature  Effects 

Friction  at  the  interfaces  between  the  deforming  preform 
and  the  dies  leads  to  a  non-homogenous  deformation  and, 
consequently,  results  in  non-uniform  densif ication.  The 
effect  of  lubrication  on  densif ication  is  shown  in  the 
macrographs,  Figure  8,  for  an  upset  forged  cylinder.  In  the 
unlubricated  case,  the  void  patterns  are  markedly  different 
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stress  is  increased  due  to  die  chill  and  it  becomes  more 
difficult  to  close  up  the  voids  with  increasing  pressure.  As 
a  result,  it  is  common  to  find  residual  porosity  at  the 
surfaces  of  hot  forged  components.  A  provision  for  good 
lubrication  is  beneficial  in  the  sense  that  it  allows  for  a 
more  uniform  densi f ication  and  pore  flattening. 

Good  lubrication  also  results  in  lower  forging  loads 
during  the  upset  forging  stage.  From  a  workability  point  of 
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Figure  8.  Macrographs  of  densif ication  characteristics  of  601AB 

aluminum  alloy  powder  preforms  upset  without  and  with 
lubrication. (3) 
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view,  lubrication  affects  the  strain  path  which,  in  turn, 
determines  the  deformation  to  fracture.  Good  lubrication 
leads  to  larger  deformation  to  fracture. 

2.7  Discussion 

The  aim  of  a  successful  powder  forging  process  is 
towards  the  production  of  a  uniformily  densified,  defect  free 
part  with  the  available  press  capacity.  A  brief  summary  of 
the  process  and  the  parameters  were  described  in  the  earlier 
sections.  It  is  evident  that  the  complexity  of  the  forging 
process,  with  its  complex  interplay  of  parameters,  precludes 
exact  quantification.  As  a  result,  considerable  experience  is 
required  for  process  designers  to  achieve  their  goal.  Proper 
preform  design  is  the  key  to  a  successful  process.  Guidelines 
and  recommendations  have  been  generated  over  the  years,  and 
these  form  the  basis  for  the  present  and  future  work. 

Prior  to  the  actual  design  of  a  preform,  prerequisite 
information  is  required  regarding  its  compaction  from  loose 
powder.  Mechanical  property  requirements  in  the  finished 
part  suggest  a  minimum  initial  preform  density.  The 
experimental  evidence,  referenced  in  section  2.2,  indicate 
that  for  low-alloy  steel  powders,  an  initial  preform  density 
of  78-80%  theoretical  is  optimum.  Compaction  of  the  preform 
to  the  required  density  may  be  achieved  either  by  die 
compaction  or  isostatic  compaction,  the  choice  of  the 
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compaction  process  being  determined  by  the  geometry  of  the 
preform.  For  die  compaction,  compaction  pressures  in  the 
range  of  40-50tsi^)  are  required.  This  requirement,  together 
with  the  available  press  capacity,  limit  the  plan  area  of  the 
preform  that  may  be  compacted. 

Prevention  of  fracture  during  powder  forging  is  the 
prime  consideration  for  preform  design.  Proper  preform  shape 
and  size  is  required  to  obtain  a  controlled  metal  flow  to 
avoid  fracture  and  undensified  regions.  In  the  case  of  hot 
repressing,  lateral  flow  is  negligible  and  fracture  is  not  a 
severe  problem.  However,  the  distribution  of  mass  in  the 
preform,  especially  in  regions  of  changing  geometry,  is 
critical  to  achieve  a  uniformily  densified  part. 

Considerable  tr ial-and-er ror  effort  is  required  to  determine 
this  mass  distribution  for  a  complex  part.^)  In  powder 
forging,  the  presence  of  flow  and  the  possibility  of  fracture 
add  considerably  to  the  uncer tainities .  Shape  and  size  of  the 
preform  become  interrelated  where  the  shape  controls  the 
fracture  likelihood  and  the  size  controls  the  mass 
distribution. 

The  workability  analysis,  in  principle,  has  been 
successfully  applied  to  preform  design  for  simple 
configurations.^)  The  material  function,  *  f ' ,  was 
determined  from  upset  tests  conducted  at  the  desired  forging 
temperature.  The  process  function,  'g',  however,  was 
experimentally  addressed  by  strain  measurements  at  potential 
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fracture  sites.  Fracture  likelihood  was  then  inferred. 

A  quantitative  prediction  of  fracture,  with  respect  to 
process  parameters,  is  a  complex  task.  A  detailed  analytical 
simulation  of  the  deformation,  together  with  a  fracture 
criterion,  at  potential  sites,  appears  to  be  the  approach. 

To  this  end,  the  state  of  analytical  development  is  reviewed 
and  assessed  in  the  following  chapter. 
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3.0  ANALYTICAL  TECHNIQUES  FOR  FORGING  OF 
SINTERED  POWDER  PREFORMS 


3.1  Introduction 


In  the  development  of  the  theoretical  mechanics  of 
plastic  deformation  of  any  material,  it  is  first  necessary  to 
establish  a  criterion  of  yielding  and  then  a  flow  rule  from 
which  the  relations  between  stress  and  strain  are  derived. 

For  plastic  deformation  of  a  conventional,  full  density,  cast 
and  wrought  material,  the  von  Mises  yield  criterion  and  the 
resulting  Levy-Mises  stress-strain  equations  are  applicable. 

Densif ication  accompanying  deformation  of  a  sintered 
powder  material  leads  to  the  violation  of  the  fundamental 
assumption  of  incompressibility  of  classical  plasticity 
theory.  Therefore,  the  yield  criterion  for  porous  materials 
must  include  the  effect  of  hydrostatic  stress. 

Kuhn  et  al  (12)  andshima*11)  derived  the  yield  criterion 
following  an  experimental  approach.  Shima  assumed  the  yield 
function  for  a  sintered  porous  material  to  be  a  linear 
combination  of  the  yield  function  of  the  incompressible 
matrix  and  the  hydrostatic  stress.  The  constants  were 
determined  from  uniaxial  compression  tests  and  correlated  to 
the  relative  density,  P  ,  of  the  porous  compact.  From 
studies  on  uniaxial  compression  tests,  Kuhn*14)  observed  a 
correlation  between  the  plastic  Poisson's  ratio,  v  ,  and  the 
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relative  density  P  .  The  empirical  relation 

v  =  0. 5Pa  (3-D 

was  exhibited  for  the  density  range  0.7<P<1.0.  The  exponent 
'a'  varied  between  1.92  and  2.0,  and  was  insensitive  to 
material  and  temperature.  As  densif ication  progresses,  the 
plastic  Poisson  ratio  approaches  0.5,  which  is  the  value  for 
conventional  fully  dense  materials.  The  yield  criterion  was 
presented  as 

f  =  {  3J^  -  (  1  -2v  )  J2  }  *  (3-2) 

A  comparitive  study  of  various  yield  criteria  published 
on  porous  materials  was  presented  by  Honess.  This  study 

demonstrated  that  equation  (3—2)  correlated  best  with  various 
experiments.  On  this  basis,  the  present  analytical  effort 
exclusively  uses  equations  (3-1)  and  (3-2). 

3.2  Basic  Equations 

3.2.1  Yield  Behavior 

The  yield  criterion  represented  by  equation  (3-2),  in 
tnrue  dimensional  stress  space,  is  an  ellipsoid  whose  major 
axis  coincides  .. i t Ii  the  direction  of  the  hydrostatic  stress 
vector.  The  intersection  of  the  yield  surface  with  the 
u-plane  is  a  circle.  As  densi f ication  progresses,  the  yield 
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surface  expands  predominantly  along  the  hydrostatic  stress 
direction.  In  the  limit,  when  v  approaches  0.5,  the 
ellipsoid  degenerates  to  the  classical  Mises  cylinder.  This 
is  schematically  illustrated  in  Figure  9. 

3.2.2  Flow  Rule 


In  conventional  plasticity  theory,  the  flow  rule  is 
obtained  from  the  normality  between  the  incremental  strain 
vector  and  the  yield  surface.  Assuming  that  this  is  valid  in 
the  case  of  sintered  powder  materials,  we  have 

deij  =  dX  .  df/d^ij  (3-3) 

Substituting  (3-2)  into  (3-3),  in  component  form,  we 

have 

d£^  =  (dX/f)  {  o^  -  v(  022  °33  )  ) 

d^22  =  (dX/f)  {  O22  —  v(  O33  +  on  )}  (3  — 2a) 

de33  =  (dX/f)  {  O33  -  v(  on  +  O22  )) 


The  term  (dX/f)  is  evaluated  by  defining  the  equivalent 
stress,  o  ,  and  the  incremental  equivalent  strain,  de  ,  for 
the  porous  compact.  Following  Hill, (1^)  effective  stress 

is  given  by 


{  U+v).  {  (o11-o22)2+(o22-o33)2+(o33-o11)?  }  +  ^2v(o11+o22+o33)2)(3-4) 
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Figure  9.  Yield  locus  predicted  by  plasticity  theory  for  porous 

materials.  Axes  represent  the  second  invariant  of  the 
stress  deviator(shear)  and  the  first  invariant  of  the 
stress  tensor (hydrostat ic  pressure).  (3) 
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The  plastic  work  done  per  unit  volume  of  the  compact  is 


(3-5) 


Using  (3-3a)  and  (3-4),  the  incremental  equivalent 
strain  is  obtained  as 


(3-6) 


and 


(3-7) 


(dA/f)  =  (de/  o) 


3.2.3  Densif ication 

The  densif ication  of  the  porous  compact  is  determined 
through  conservation  of  mass.  In  general,  during 
deformation,  density  gradients  may  develop  in  the  domain  of 
the  compact.  The  equation  of  conservation  of  mass  is  given 
by 


3p/3t  +  3p/  ^  +  p3uk/3xk  =  0 

where  uk  :  components  of  the  velocity  vector 
xk :coordinates  of  a  point. 


(3-0) 


k  :is  the  summation  index. 
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3.3  Application  to  Simple  deformations 

The  plasticity  theory  developed  above  has  been  applied 
to  predict  the  densi f ication  and  deformation  of  a  porous 
compact  in  simple  applications,  namely,  frictionless  uniaxial 
compression,  plane  strain  compression,  repressing  and 
hydrostatic  compression.  The  results  are  presented  in 
Reference  3,  and  good  predictive  capability  was  demonstrated, 
establishing  the  soundness  of  the  theory  and  its  underlying 
principles.  In  the  above  applications,  frictional  constraint 


at  the 

die  surfaces 

is 

zero  or 

does  not 

influence 

the 

deforma 

ition  hence, 

the 

densif i 

cation  is 

uniform. 

As 

result. 

equation  (3 

-8) 

may  be 

simplified 

and  cast 

in 

incremental  form 


-dp/p  =  de^  +  de22  +  df:33  (3-8a) 

The  govering  equations  can  be  integrated  in  closed  form 
with  this  simplification.  However,  when  friction  at  the  die 
surfaces  contributes  to  non-uniform  densif ication ,  a  closed 
form  solution  to  simple  deformations  is  not  feasible. 
Numerical  techniques  have  to  be  resorted  to,  and  this 
constitutes  the  main  effort  of  this  chapter. 


3.4  Approximate  Techniques 


Exact  solutions  to  conventional  plasticity  problems  are 
lacking  except  for  very  simple,  idealized  deformations. 
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Approximate  techniques  were  developed,  which  are  listed  in 
order  of  increasing  sophistication  and  effort,  as  follows: 

(1)  Slab  or  Equilibrium  Analysis. 

(2)  Upper  Bound  Techniques. 

(3)  Finite  Element  Method. 

The  slab  analysis  can  predict  the  effect  of  interface 
friction  and  other  process  parameters  on  the  total  forging 
load.  The  upper  bound  technique  minimizes  the  deformation 
energy  resulting  from  an  assumed,  kinematically  admissible, 
velocity  field  and  yields  a  closer  approximation  to  the  total 
forging  load.  Both  techniques  are  incapable  of  detailed 
prediction  of  stress  and  strain  distributions.  The  finite 
element  method  can  predict  distributions  of  stress  and  strain 
in  the  deforming  material.  However,  applications  are 
presently  limited  to  simple  geometries  due  to  the  extensive 
computational  effort  and  complexities  in  the  numerical 
procedures.  This  method  is  currently  being  researched 
intensely . 

With  the  objective  of  developing  comparable  approximate 
techniques  for  sintered  powder  materials,  two  techniques  are 
presented.  The  first  is  the  slab  analysis,  which  is  capable 
of  predicting  approximations  to  the  total  load  and  density 
distributions  in  the  presence  of  friction  at  the  die  walls. 
The  second  technique  is  a  hybrid  scheme,  based  on  an  energy 
bounding  principle  but,  numerically  similar  to  the  finite 
element  method.  This  scheme  was  first  introduced,  for 
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conventional  fully  dense  materials,  by  Lee  and  Kobayashi . ^ 
Applications  are  presently  limited  to  simple  geometries  for 
purposes  of  evaluation. 

3.5  Slab  Analysis  for  Axisymetric  Upset  Forging 
3.5.1  Development  of  Equations 

The  underlying  assumptions  for  the  slab  analysis 
are:*17* 

(i)  Stress  on  a  plane  or  surface  normal  to  the  flow 
direction  is  a  principal  stress.  Hence,  for  an  infinitesmal 
element  parallel  to  this  plane,  the  deformation  element  is 
considered  homogenous. 

(ii)  The  presence  of  external  friction  or  shear  does  not 
modify  the  internal  stress  distribution. 

Consider  the  upsetting  of  a  compact  whose  cross-section, 
shown  in  Figure  10,  is  a  trapezoid.  For  an  infinitesmal 
element,  the  radial  stress,  °r  ,  is  developed  due  to  the 
frictional  shear  on  the  punch  interfaces.  Under  the 
assumption  of  homogenous  deformation,  °r  ,  does  not  vary  with 
axial  position.  Assuming  that  the  hoop  stress,  a0  ,  is  equal 
to  the  radial  stress,  the  yield  criterion,  (3-2),  simplifies 
to 

Y 2  =  (o  -  a  )2  +  (l-2v)  (a  +  2a  a  )  .  (3-9) 

rz  r  r  z  '  * 

where  Y  =  yield  strength  of  the  compact  at  density  P  . 
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(A)  Trapezoidal  cross  section  showing  zone  subdivision. 

(B)  Geometry  of  nth  zone  with  infinitesmal  element. 

(C)  Stresses  acting  on  infinitesmal  element. 


Figure  10.  Geometry  of  axisymmetric  component  with  a  trapezoidal 
cross-section  for  Slab  analysis. 

I 
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°2  =  applied  normal  stress. 

The  flow  equations,  (3-3a) ,  expressed  in  terms  of  strain 
rates  are: 


e  =  (  \/  Y)  (o  -  v (o  +  a  )  } 
r  r  r  z 

(3-10) 

e  =  (  X/  Y)  {a  -  2 vo  } 
z  z  r 

The  radial  velocity,  ,  at  the  radial  position  r,  is 


ur  =  uo  +  f  K  {  K/K)dr  (3-ID 

ro 

where  uQ=  radial  velocity  at  r  =r0. 

The  densif ication  rate,  (3-8),  assuming  density  as  a 
function  of  radial  position,  is 


8p/9t  =-  pe{l+2(e/e)}+  u  9p/3r 

z  r  z  r 


(3-12) 


The  distribution  of  the  radial  stress,  ar ,  is  obtained 
from  equilibrium  considerations  of  an  infinitesmal  element. 
From  Appendix  A,  we  have 


do  /dr  =  -{(a  -a  ) (tana  +tana,)  -  (x  Sec2a  +x,Sec2a,)  }  /  h 
r  zr  u  1  u  ul  1 


=  r  ;  a  =  (a  ) 
n+1  r  r  r=r 


n+1 


r 


(3-13) 
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3.5.2  Solution  of  Equations 

Equations  (3-9)  to  (3-13)  together  with  (3-1)  form  the 

set  of  equations  for  an  incremental  numerical  procedure. 

This  was  implemented  on  a  digital  computer  and  the  upset  of  a 

circular  cylinder  was  simulated  for  various  values  of  aspect 

ratios  and  interface  friction. 

The  cross-section  was  descretized  into  n  slabs.  Density 

and  velocities  are  specified  at  the  sides  flanking  each  slab. 

The  density  is  assumed  to  vary  linearly  within  each  element. 

For  a  given  axial  strain-rate,  ,  assumed  constant,  equation 

(3-13)  coupled  with  (3-9)  is  numerically  integrated  from  the 

outermost  slab  towards  the  interior  by  a  fourth  order 

Runge-Kutta  scheme.  Eqaution  (3-9)  then  furnishes  a  ,  whence 

z 

the  ratio  ( e  /  H  is  evaluated  from  (3-10).  Integrating 
(3-11),  the  radial  velocity  distribution  for  the  current 
increment  is  obtained  with  the  densif ication  computed  from 
(3-12).  For  a  time  increment  At,  the  density  and  positions 
of  the  slabs  are  updated  and  the  procedure  repeated. 

Numerical  instabilities  were  encountered  in  the  above 
procedure  for  the  cases  of  high  interface  and/or  low  aspect 
ratios.  For  low  values  of  friction,  the  procedure  was  stable 
and  yielded  results  consistent  with  experimental  trends.  This 
prompted  a  closer  look  on  the  influence  of  friction  on  the 
algorithm. 
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For  a  circular  cylinder,  au  =  a1 =  0.  Assuming  that  Tu  =T 
=mk,  where  m  is  the  friction  factor  and  k  the  shear  yield 
strength  of  the  compact,  equation  (3-13)  may  be  integrated  to 
yield 


a  /  Y  =  -{  m/ (H/D) }  {1  -  (r/R) }  /  /3 
r 


(3-14) 


It  is  evident  that  (ar/  Y  )  is  directly  proportional  to 
'm'  and  inversely  to  (H/D).  High  values  of  m  and  low  values 
of  (H/D)  result  in  large  values  of  (  °r/  Y),  these  being  the 
situations  when  the  algorithm  is  unstable. 

Solving  the  yield  criterion,  (3-9),  explicitly  for 
(  0  /  Y ) ,  we  have , 

o  /  Y  =  2vor/  Y  -  {  1  -2(ar/  Y) 2 (1+v)  (l-2v)  }*  (3-15) 

This  is  plotted  in  Figure  11  for  various  values  of  v  . 

We  infer  that  as  ( °r/  Y  )  decreases  linearly  from  zero,  as 
suggested  by  (3-14),  (  a?/  Y  )  attains  a  maximum  negative 

value.  The  maximum  value  is  easily  derived  to  be 


(3-16) 


It  can  be  shown  that  (3-16)  corresponds  to  the  stresses 
cn  a  circular  disk  under  a  repress  deformation. 

The  strain-rate  ratio,  (3-10),  changes  sign  in  the 
vicinity  of  the  stress  state  (3-16)  according  to 


AXIAL  STRESS 


39 


RADIAL  STRESS  -<T"r/Y 


Figure  11.  Variation  of  radial  stress  vs.  axial  stress  for  disk  upset 
as  predicted  by  the  yield  criterion. 


i^o 


a  /  Y  <  (a  /  Y) 
r  r 

a r/  Y  (a  /  Y) 

a  /  Y  >  (a  /  Y) 
r  r 


(3-17) 


Suppose  that  the  disk  has  a  low  aspect  ratio  and  is 

deformed  under  a  large  frictional  constraint.  For  material 

flow  directed  radially  outward,  the  radial  stress  varies 

linearly  and  monotonical ly  from  zero  at  the  rim  to  a  finite 

value  at  the  center.  The  strain-rate  ratio  is  negative  at 

the  rim,  consistent  with  an  outward  flow,  and  increases 

toward  the  center.  Suppose  at  some  interior  point,  r=r *, 

( e  /c  )  is  zero.  Then  for  any  position  r<r*,(e/e  )  is 

r  z  r  z 

positive  implying  that  the  flow  is  radially  inward.  This 
would  imply  a  reversal  in  the  direction  of  interfacial  shear 
and  result  in  a  decrease  in  (-  °r//  Y  )  .  This  contradicts  the 
hypothesis  that  (-  Y)  is  always  increasing.  Hence,  we 
conclude  that  (- ar/  Y  ) i (- ar/  Y ) * ,  which  implies  that  the 
radial  stress  distribution  (3-14)  is  strictly  given  by 


a  /  Y 
r 


(m/(H/D)  }  { 1-r/Ft}  /  /3 

*  * 

(a  /  Y)  ;  0  <  r  <  r 


(3-14a) 


Physically,  this  means  that  for  extreme  case  of  porous 
parameters,  there  are  two  zones  of  deformation.  An  outer 
zone  undergoing  upset  deformation  surrounding  an  inner  zone 
undergoing  a  pure  repress  deformation.  It  can  be  shown  that 
as  the  deformation  progresses,  the  ensuing  densif ication  of 
the  repress  core  results  in  the  core  shrinking  to  zero, 
whence  the  entire  disk  then  undergoes  an  upset  deformation. 

The  inclusion  of  equation  (3-14a)  into  the  numerical 
procedure  resulted  in  a  stable  procedure  for  all  conditions 
of  process  parameters. 

3.5.2  Results  for  Disk  Upset 


Available  experimental  data  on  upsetting  of  disks  under 
different  conditions  of  lubrication  and  aspect  ratios  are 
shown  in  Figures  12  and  13.  The  overall  densif ication, 

Figure  12,  is  enhanced  with  deformation  under  high  frictional 
constraint  and  for  low  aspect  ratios,  respectively.  The 
density  gradients  at  various  levels  of  deformation,  Figure 
13,  are  consistent  with  the  macrographs  of  Figure  8. 

Results  of  simulation  by  the  slab  analysis,  under 
conditions  approximate  to  the  experiments,  are  shown  in 
Figures  14,  15  and  16.  The  friction  condition,  qualitatively 
defined  in  the  experimental  data  as  'lubricated'  and 
'unlubricated',  is  specified  in  the  analysis  by  a  constant 
interfacial  shear.  The  friction  factor,  m,  for  the  shear 
specification  was  chosen  as  (3.2  and  0.5  for  the  lubricated 
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Figure  12.  Experimental  data  on  overall  dens i f i ca t i on  in  disk  upset 
of  601AB  aluminum  alloy  powder  preforms  for  various 
aspect  ratios  and  lubrication.  (3) 
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Figure  13.  Dens i f i cat  ion  profile  as  measured  across  an  upset  60IAB 
aluminum  alloy  powder  disk  preform. (3) 
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Figure  14.  Overall  dens i f i ca t i on  calculated  by  Slab  analysis  for 

upset  of  601AB  aluminum  alloy  powder  disk  preforms  at 
various  aspect  ratios  and  lubrication. 
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Figure  15.  Comparison  of  overall  dens i f i cat i on  calculated  by  Slab 
analysis  with  experimental  data  for  upset  of  601AB 
aluminum  alloy  powder  disk  preforms. 
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Figure  16.  Comparison  of  dens i f i ca t i on  profiles  calculated  by  Slab 
analysis  with  experimental  profiles  for  upset  of  601AB 
aluminum  alloy  powder  disk  preforms. 


and  unlubricated  conditions,  respectively. 

Trends  in  the  overall  densif ication  with  deformation  as 


predicted  by  the  slab  analysis,  Figure  14,  are  consistent 
with  experimental  observation.  Figure  15  illustrates  a 
typical  comparison  with  experimental  data.  The  analysis 
underestimates  the  actual  overall  densif ication.  This  may  be 
the  result  of  an  underestimation  of  the  friction  factor. 

Figure  16  illustrates  the  comparison  of  the  density 
gradient  predicted  by  the  analysis  with  the  experimental 
data.  There  exists  a  marked  difference,  which  immediately 
raises  questions  on  the  validity  of  the  slab  analysis.  The 
data  indicates  that  the  density  is  fairly  uniform  in  the 
central  regions  and  falls  off  towards  the  rim,  whereas  the 
analysis  yields  an  almost  linear  density  distribution.  The 
uniform  density  at  the  central  regions  indicate  the  presence 
of  a  core  under  repress  deformation,  but  for  the  process 
conditions  selected  this  is  false. 

Rooyen  and  Backofen^®)  attempted  to  experimentally 
determine  the  distribution  of  interfacial  shear  for  the 
upsetting  of  Aluminum  disks.  Their  observation  was  that  the 
interfacial  shear  was  not  constant  but  increased  toward  the 
rim  of  the  disk.  The  variation  was  complex  and  not 
quantifiable. 

Variation  in  the  frictional  constraint  at  the  interface 
influences  the  densif ication  behavior  of  a  porous  compact. 
Interfacial  friction  generates  the  radial  pressure  (cf 


equation  3-14),  the  variation  of  which  results  in  density 
gradients.  If  we  assume  that  the  frictional  shear  is 
distributed  linearly  and  increasing  towards  the  periphery  of 
the  disk,  the  radial  pressure  distribution  becomes  parabolic 
(as  opposed  to  linear  for  a  constant  shear).  This  would 
result  in  the  dens  if ication  being  dampened  near  the  core  of 
the  disk,  as  suggested  by  the  data.  To  determine  the  extent 
of  the  change  that  would  result,  a  linear  shear  specification 
was  imposed.  The  shear  was  zero  at  the  center  and  with  a 
slope  such  that  the  average  shear  was  equal  in  magnitude  to 
the  constant  shear  that  was  specified  earlier. 

The  results  are  shown  in  Figure  16  by  the  broken  line. 
The  density  distribution  predicted  resembles  the  data 
closely,  especially  at  the  higher  reductions.  The  effect  on 
the  overall  dens i f ication  ,  Figure  15,  is  marginal,  but  tends 
towards  a  better  agreement  with  the  data. 

3.6  Energy  Method  for  Axisymmetric  Upset 

3.6.1  Derivation  of  the  Functional 

Consider  a  deforming  porous  body  in  a  fully  plastic 
state.  Let  denote  the  surface  where  velocities  are 
prescribed  and  Sp  the  surface  where  tractions  are  specified. 
For  this  body  of  volume  V,  let  the  following  be  defined: 


U9 


u  :  proscribed  velocity  vector  on  . 

F  :  prescribed  traction  vector  on  SF  . 

:  actual  strain-rate  field  satisfying  compatibility. 
:  equilibrium  stress  field  satisfying  boundary 
conditions. 

u*  :  any  velocity  field  that  satisfies  the  prcscribeo 
boundary  conditions. 


e* 

ij 


strain-rate  field  derived  from  u* 

i 


:  stress  field  derived  from  via  the  flow  rule 

and  satisfying  the  yield  criterion. 


eg 


eq 


n . 
D 


equivalent  stress  derived  from  components  of  . 

equivalent  strain-rate  field  derived  from  components 
of  ej  . 


components  of  unit  outward  normal  to  surface. 

The  principle  of  maximum  plastic  work  for  a  concave  -to- 
origin  yield  locus  is  given  by^^ 


/  (o* .  -  a.  .)  e* .  dv  >0 

V  ^  ^ 


(3-18) 


By  the  principle  of  virtual  work,  we  have 


/  o.  .  e*.  dV  =  f  E.u*dS  + 
V  «  sF  ' 

Using  equation  (3-18)  we  have. 


/  (a . ,n . )U. dS 
c  ID  D  i 


4>  =  fa* .  e*  .dV  -  f 

V  13  s 


F.ues  >  f  (a. . n . ) U . dS 
e  ID  D  1 
F  U 


(3-19) 
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Following  Hill^®^ 


fa  e. .dv  =  /  F.udS  +  /  (a..n.)U.  ds 

V  13  13  S*'  S  13  3  1 

F  u 


(3-20) 


Hence,  from  (3-19)  and  (3-20) 


I  o * . e * ,  dV  -  /  F.u*  dS  >  /  a, .e.  ,  dv  - 

»  ««  s''  V  13  13 

F 

=  /  (o  .n  )U .  dS 

s  13  j  1 

U 


5  F.u  dS 
S„  ' 


The  functional,  in  terms  of  equivalent  stress  and 
strain-rate,  is  then  given  by 


4>  „  /  a  e  dV  -  /  F .  u  ds 

v  eq  eq  Sp  ' 

3.6.2  Solution  of  Equations 

We  seek  to  determine  the  kinematically  admissible 
velocity  field  u*  ,  from  a  class  of  velocity  fields 
satisfying  the  prescribed  boundary  conditions,  that  minimizes 
the  functional  <t>  .  The  stress  field  a*..  ,  then  is  an  upper 
bound  on  the  actual  stress  field  o^,  as  o*_,  need  not 
necessarily  satisfy  the  equilibrium  conditions. 

To  obtain  a  solution  for  the  kinematically  admissible 
velocity  that  minimizes  the  functional  *  ,  we  use  the  finite 

element  procedure  of  Lee  and  Kobayashi.  The 

axisymmetric  domain  is  descretized  into  M"  elements.  The 
velocity  field,  u* ,  is  approximated  in  terms  of  the 
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velocities  at  the  nodal  points,  y,  by 


u*  =  Gv 


(3-22) 


The  strain-rate  field  is  derived  as 


e 


=  Bu* 


(3-23) 


The  functional  <t>  is  approximated  as 


m  ,  M 
=  E  4, =E 


m=l 


{  /  o  e  dV 
m=l  v(»)^  eq 


-  f  F.u*  dS  } 
„  (m) 


(3-24) 


The  equivalent  strain-rate,  eeq  ,  is  then  derivable  as 


•  o  •  T  • 

cz  =  2  e  Ae  + 

6q  3(l+v) 


•  x  • 

e  Ie 


3(l-2v) 


(3-25) 


where  I  is  the  identity  matrix,  A  is  a  constant  matrix. 
Substituting  (3-22,  23,  24)  into  (3-25),  we  have 


M  i 

4>  =  I :  {  /  O  (yTXv)J  dV  -  /  FTGy  dS  }  (3-26) 

m=l  (m)  eq  (m) 

V  SF 

where  x  =  btab  +  btib 

Minimizing  4  with  respect  to  the  nodal  point  velocities. 


M  , 

=  Z  {  /  o  (yTXv)-4  Xv  dV  -  f  GT£  dS  } 
m=l  (m)  ^  (m) 

V  SF 


3^/dV 


(3-27) 
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This  represents  a  system  of  non-linear  equations  in  the 
nodal  point  velocities.  The  solution  determines  the  velocity 
field  that  minimizes  the  functional.  This  solution  may  be 
obtained  either  by  the  Newton-Raphson  method  or  by 
linearization  of  equation  (3-27). The  latter  is 
accomplished  by  assuming  that  at  the  n-th  iteration,  the 
velocity  vector  is  given  by 


v 


(n) 


Y(n-1) 


+  Ay 


(n) 


(3-28) 


The  linearization  procedure  and  the  resulting  linear 
system  of  equations  in  Ay(  is  detailed  in  Appendix  B. 
Schematically,  the  system  of  equations  is 


M 


l 

m=l 


a  P 
eq 


(n-1) 


Av,  , 
-  (n) 


M 

E 

m=l 


a  H, 
eq- (n-1) 


M 

-IT 

m--l 


=  0 


(3-29) 


The  iteration  procedure  for  an  increment  of  deformation 
is  started  by  assuming  a  trial  velocity  field.  This  is 
iterated  till  the  velocity  iterate,  Av  t  is  less  than  a 
prescribed  tolerance. 

Once  a  converged  velocity  field  is  obtained  for  an 
increment,  the  strain-rate  distribution  may  be  obtained  from 
(3-23).  Stress  distribution  is  obtained  by  inverting  the 
flow  rule.  Note,  however,  that  this  inversion  is  valid  only 
for  sintered  powder  materials. 
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3.6.3  Results  for  Disk  Upset 

As  an  initial  application  of  the  energy  method,  the  disk 
is  descretized  into  elements  which  are  slabs.  This  element 
configuration  was  chosen  to  provide  a  check  for  the  results 
of  the  slab  analysis.  For  a  more  detailed  analysis,  the 
element  may  be  a  quadrilateral,^*^  in  which  case  even  the 
attendent  bulging  of  the  free  surface  may  be  obtained. 

The  element  equations  and  material  characterization  are 
given  in  Appendix  B.  The  velocity  field  is  approximated  by  a 
linear  interpolation  of  the  velocities  of  the  element  sides. 
Density  is  also  assumed  to  vary  linearly  within  the  element. 
For  the  first  deformation  increment,  the  assumed  trial 
velocity  field  was  that  of  uniaxial  frictionless  compression. 
For  subsequent  increments,  the  converged  velocity  field  of 
the  previous  increments  is  used.  The  number  of  iterations 
required  to  meet  the  specified  error  tolerance  depend  upon 
how  close  the  assumed  velocity  field  is  to  the  actual 
velocity  field.  The  convergence  criteria  was  specified  as 

I  I A  Y I |/| | v | I  <  1.0X10"4 

where  IIavII,  I  I y I  I  are  the  Euclidean  vector  norms.  On  the 
average,  a  converged  velocity  field  was  obtained  in  6-8 
iterations  for  the  first  increment,  and  the  subsequent 
iterations  required  2-3  iterations. 

The  deformation  of  the  disk  was  simulated  under  the  same 
conditions  that  were  used  in  the  slab  analysis.  Interface 


shear  was  specified  as  constant  in  one  set  of  computations, 
and  varying  linearly  in  another  set.  The  computed  density 
distributions  for  both  shear  specifications  are  shown  in 
Figure  17  along  with  the  experimental  data. 

The  predictions  are  similar,  both  qualitatively  as  well 
as  quantitatively,  to  that  of  the  slab  analysis.  This  result 
is  not  surprising  for  the  type  of  element  chosen.  For  the 
case  of  zero  interfacial  shear,  both  techniques  yield  results 
that  are  identical  and  agree  with  the  results  for  uniaxial 
frictionless  compression. 

The  predicted  overall  densif ication  is  also  similar. 

The  predicted  forging  load  is  slightly  higher  for  the  energy 
method  when  compared  to  that  of  the  slab  method.  The  good 
correlation  between  the  results  for  the  same  set  of 
conditions  establishes  the  correctness  of  the  algorithms  and 
tne  computations. 


3.7  Discussion 

Two  numerical  techniques,  applicable  to  sintered  powder 
materials,  were  presented  in  the  above  sections  and  applied 
to  the  uniaxial  compression  of  a  cylinder  with  interfacial 
frictional  constraint.  The  predicted  results  are  consistent 
with  experimental  observations.  The  pros  and  cons  of  both 
techniques  are  compared  with  respect  to  accuracy, 
computational  effort  and  limitations.  The  candidacy  of  both 
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NORMALIZED  POSITION  r/R 


Figure  17.  Comparison  of  dens i f i cat  ion  profiles  calculated  by  Energy 
method  with  experimental  profiles  for  upset  of  601AB 
aluminum  alloy  powder  disk  preforms. 
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for  analytical  process  modeling  is  discussed  in  the  following 

f aragraphs  . 

The  slab  analysis,  by  virtue  of  its  underlying 
assumptions,  is  incapable  of  predicting  the  local  states  of 
stress  and  strain  during  deformation.  In  the  case  of  wrought 
materials,  the  primary  advantage  is  in  obtaining  quick 
estimates  of  the  forging  load.  The  assumption  of  homogenous 
deformation,  valid  strictly  in  the  case  of  frictionless 
compression,  is  not  valid  in  reality.  This  assumption, 
however,  considerably  simplifies  the  mathematics  required  to 
obtain  an  expression  for  the  forging  load.  In  the  case  of 
sintered  powder  materials,  the  resulting  non-uniform 
densif ication  precludes  the  possibility  of  an  explicit 
expression  for  the  forging  load. 

The  energy  method,  on  the  other  hand,  does  not  assume  a 
homogenous  deformation.  Yet,  the  results  of  the  slab  analysis 
and  the  energy  method  compare  very  well.  Figure  18 
illustrates  the  typical  computed  radial  and  hoop  strain-rates 
for  the  slab  analysis  and  the  energy  method.  The  assumption 
of  homogenous  deformation  seems  to  lead  to  strain-rates  that 
are  an  average  of  the  radial  and  hoop  strain-rate  values  of 
the  energy  method. 

Ihe  energy  method  has  the  potential  ability  to  describe 
the  stress  and  strain  distribution  within  a  deforming 
material.  The  present  element  chosen,  however,  results  in  a 
considerable  oversimplification  and  as  a  result  does 
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NORMALIZED  POSITION  r/R 


Figure  18.  Comparison  of  radial  and  hoop  strain-rates  calculated  from 
the  Energy  method  with  strain-rates  of  the  Slab  analysis. 
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injustice  to  the  capability  of  the  technique.  As  mentioned 
earlier,  this  element  was  chosen  primarily  to  evaluate  the 
merit  of  the  technique  and  provide  for  evaluating  the  slab 
analysis.  A  quadrilateral  element,  with  bilinear  velocity 
interpolation,  is  the  simplest  configuration  that  would 
provide  the  detailed  description  of  the  deformation.  This 
extension  is  a  subject  for  future  activity. 

To  compare  the  computaional  effort  and  accuracy  of  the 
techniques,  the  process  conditions  were  made  similar.  Ten 
slabs  were  used  in  both  cases  and  the  increment  of 
deformation  was  chosen  to  be  equivalent  to  an  axial  strain  of 
0.01.  This  choice  is  a  trade  off  between  the  computational 
effort  and  errors,  to  reach  a  specific  height  reduction. 

For  the  slab  analysis,  errors  were  measured  by  the 
difference  between  the  computed  mass  at  each  increment  and 
the  original  mass  of  the  cylinderical  preform.  For  the 
deformation  increment  chosen,  the  error  was  held  to  less  than 
1%  ,  for  a  height  reduction  of  50%.  Increasing  the  number  of 
slabs  does  not  contribute  in  any  positive  way  to  the  results, 
but  only  increases  the  computational  effort.  Hence,  some 
judgement  needs  to  be  exercised  regarding  the  choice  of  the 
deformation  increment  and  the  number  of  slabs.  The  slab 
analysis  is  limited  in  its  capability  of  process  modeling. 

The  energy  method,  being  an  iterative  process,  is  by  far 
a  better  technique.  For  the  slab  type  element  chosen,  the 
computational  effort  is  approximately  of  the  same  order.  No 
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attempt  was  made  to  optimize  the  number  of  elements  as  it  is 
recognized  that  the  slab  element  does  not  exploit  the 
capability  of  this  method.  It  needs  to  be  cautioned  that  in 
going  from  a  slab  type  element  to  a  quadrilateral  element 
will  resul't  in  a  large  increase  in  computation.  This 
increase  can  be  justified  by  the  detailed  results  that  is 
potentially  obtainable. 

When  full  densif ication  is  approached,  numerical 
difficulties  may  arise.  Specifically,  this  occurs  because  as 
full  density  is  approached,  the  Poisson's  ratio  tends  to  0.5 
and  some  terms  in  the  equations  may  become  very  large.  This 
problem  can  be  avoided  if  full  density  is  assumed  to  be 
99.95%  theoretical.  It  was  numerically  verified  that  as  full 
density  is  approached,  the  volumetric  strain-rate  tends  to 
zero,  i.  e.,  the  incompressibility  condition  is  approximately 
satisfied.  The  significance  of  this  is  clear  when  one 
compares  the  current  method  for  sintered  powder  materials  at 
99.95%  density  with  the  companion  method  for  wrought 
materials  of  Lee  and  Kobayashi.  For  wrought  materials, 

the  incompressibility  condition  is  absorbed  into  the 
functional  by  means  of  a  Lagrange  multiplier.  This  is 
computed  for  each  element,  for  each  iteration,  and  increases 
the  number  of  equations  to  be  solved  by  the  number  of 
elements  in  the  descretized  domain.  The  Lagrange  multiplier 
was  shown  to  be  equivalent  to  the  hydrostatic  stress.  To 
compute  the  stress  distribution,  the  deviatoric  stresses  are 
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calculated  from  the  strain-rates  and  the  hydrostatic  stress 
is  added.  For  sintered  powder  materials,  the  stress 
distribution  is  obtainable  by  directly  inverting  the  flow 
rule.  Thus,  for  large  scale  problems,  considerable  savings 
in  computations  and  as  well  as  storage  can  be  realized.  The 
accuracy  deficiencies,  if  any,  needs  to  be  assessed. 

An  interesting  result  of  the  present  analytical  effort 
was  the  influence  of  interfacial  friction  on  the  process 
modeling.  Considerable  research  was  done  in  the  last  two 
decades  to  understand  and  quantify  the  effect  of  friction  in 
metalworking  processes.  Current  practice  is  to  assume  the 
loti  'odd  or  the  constant  friction  factor  (Amoton's  Law), 
l  latter  being  preferred  because  of  the  mathematical 
simplifications  that  result  without  appreciable  differences 
in  the  forging  load.  A  quantitative  average  estimate  of  the 
friction  factor  is  obtained  by  the  ring  test,^l)  where 
changes  in  the  geometry  of  the  ring  are  used  as  an  indicator 
of  the  magnitude  of  friction.  However,  it  was  noted  in  the 
present  work  that  changing  the  frictional  shear  distribution 
has  a  pronounced  influence  on  the  densif ication  profiles  but 
virtually  no  change  in  the  deformed  geometry.  The  change  in 
the  forging  load  was  also  marginal,  as  shown  in  Figure  19. 
Perhaps  the  ring  test  for  a  sintered  powder  material  would 
shed  more  light  on  a  technique  for  assessing  friction 
factors.  The  energy  method,  with  the  slab  type  element  would 
be  a  good  approach  to  take  to  assess  ring  deformation. 


LOAD  ( I  hs ) / RAD  I  AN 
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7  HEIGHT  REDUCTION 


Figure  19.  Calculated  forging  load  for  upset  of  60IAB  aluminum  alloy 
powder  disk  preforms. 
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4.0  COMPUTER-AIDED  DESIGN  OF  PREFORMS 

4.1  Introduction 

The  importance  of  preform  design  was  realized  at  a  very 
early  stage  in  the  development  of  powder  forging.  An  early 
example  is  the  forging  of  differential  side  pinion  gears, ^2) 
where  a  partial  bevel  preform,  determined  from  experimental 
results  using  plasticine,  successfully  led  to  a  sound 
forging.  Subsequently,  the  influence  of  preform  shape  on 
material  flow,  residual  porosity  and  occurrence  of  flaws  was 
further  emphasized  in  the  forging  of  specific  axisymmetric 

/on  O  1  \ 

parts.'  J •  q)  An  experimental  approach  was  pursued  in  these 
works,  where  the  actual  preform  was  developed  after  a  study 
of  the  mechanics  of  deformation  and  defect  formation  in  a 
limited  series  of  model  experiments.  Kuhn  et.  al . ^ •  25-27) 
pursued  this  approach  in  some  detail,  the  result  of  which  is 
the  establishment  of  preform  design  as  the  single  most 
important  parameter  in  powder  forging.  Empirical  guidelines 
for  the  prevention  of  defects  for  a  variety  of  simple  shapes 
were  also  obtained. 

The  application  of  these  concepts  to  a  complex  forging 
has  been  successful,  but  is  quite  tedious. ^  Firstly,  an 
understanding  of  metal  flow,  densif ication  and  fracture 
behavior  of  preforms,  acquired  primarily  with  experience,  is 
required  to  determine  the  shape  of  the  preform.  Secondly, 
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the  preform  is  to  be  sized  so  as  to  maintain  the  desired 
shape ,  provide  adequate  mass  distribution  in  various  regions 
and  maintain  the  proper  weight  control.  The  latter  task  is 
time  consuming  and  prone  to  human  errors,  particularly  in  the 
case  of  complex  shapes. 

Accurate  determination  of  volumes  lends  itself  easily  to 
computerization.  The  computerization  of  the  mechanics, 
however,  is  not  straight  forward.  If  this  is  achieved,  the 
oesign  of  powder  preform  forging  processes  may  be 
accomplished  by  those  who  are  unfamiliar  with  the  details. 
Development  of  the  computer-aided  design  approach  was 
motivated  by  this  need.  This  original  contribution  is  aimed 
towards  exploitation  of  powder  forging  as  the  technology  for 
the  new  decade. 

4.2  Review  of  CAD  in  Manufacturing 

The  use  of  computers  in  manufacturing  is  growing 
rapidly,  primarily  due  to  increased  productivity  and  superior 
quality  control  that  is  achieved.  The  earliest  use  of 
computers  was  in  Numerical  Control (NC)  Machining.  The  NC 
language,  APT  (Automated  Programmed  Tool),  was  developed  for 
this  purpose.  APT  is  used  to  define  the  part  geometry  and  the 
necessary  machine  tool  control  commands  to  drive  a  cutting 
tool  to  produce  the  part.  APT  developed  over  a  decade  to 
encompass  more  complex  part  geometries;  requiring  trained 
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personnel  to  implement  its  sophisticated  capabilities.  In 
parallel,  machine  tool  builders  developed  the  necessary 
hardware  to  complement  the  advances.  APT  was  designed  to  aid 
machining  technology,  and  it  ably  does  so. 

The  application  of  similar  technology  to  precision 
forging  was  attempted  by  researchers  at  Battelle  Columbus 
Laboratories.^^  a  modular  approach  is  followed,  where 

representative  cross-sections  are  isolated;  preforms  are 
designed  using  established  practices  and  then  assembled  in  a 
building  block  fashion  to  yield  the  preform  for  a  complex 
part.  One  approach  is  limited  to  a  class  of  structural 
components  of  <  rib-web  type,  where  any  representative 
cross-section  is  an  assembly  of  basic  L  shapes.  The  preform 
for  a  L  shape  is  defined  as  having  a  specific  but  different  L 
shape,  which  is  sized  to  maintain  a  constant  area  equal  to 
the  area  of  the  finished  L.  Limited  or  no  user  intervention 
is  required  in  this  design  scheme.  The  analysis  is  done  on  a 
cross-section  basis,  primarily  because  analytical  techniques 
are  currently  only  available  for  plane-strain  or  axisymmetric 
deformations.  The  major  obstacle,  in  the  work,  is  the  part 
description  and  the  determination  of  representative 
cross— sections. _  APT  is  used  for  this  purpose. 

Preform  design  for  powder  forging  has  a  different  set  of 
rules  compared  to  those  for  conventional  forging.  For 
non-axisymmetr ic  structural  components,  the  use  of  the  L 
shaped  mociule  is  not  feasible  due  to  unequal  densif ication 


66 


that  occurs  in  the  rib  and  web.  The  current  tr ial-and-er ror 
practice  of  preform  design  identifies  "regions"  in  the  part, 
and  the  preform  is  obtained  by  distributing  the  material  in 
various  regions.  The  density  of  the  material  in  the  regions 
need  not  be  the  same,  depending  on  whether  the  forging  is 
accomplished  by  repressing  or  not. 

The  earlier  chapters  summarized  the  powder  forging 
process  and  the  parameters  that  affect  the  process.  The 
state  of  analytical  techniques  was  also  examined,  for 
possible  application  to  the  current  work.  In  the  following 
sections,  the  question  of  computer-aided  preform  design  for 
powder  forging  is  addressed. 

4.3  Basic  Considerations  for  Preform  Design 

Structural  forgings,  in  general,  have  complex  geometries 
which  result  in  complex  metal  flow.  Certain  elementary  or 
generic  metal  flow  patterns  can  be  identified,  which 
collectively  result  in  a  complex  flow.  Each  of  these  generic 
flow  patterns  are  result  of  a  unique  combination  of  the 
preform  shape  with  respect  to  the  die  profiles. 

Each  generic  flow  pattern  has  associated  with  it  unique 
characteristics  regarding  defect  and  fracture  likelihood, 
ease  of  forging  etc.  Once  these  generic  elements  are 
established  in  a  complex  forging,  a  rational  preform  design 
procedure  may  be  effected. 
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4.3.1  Forging  Direction  Specification 

The  choice  of  forging  direction  affects  the 
producibility  of  the  part  by  powder  forging.  In  general,  the 
choice  must  be  such  that  the  part  may  be  ejected  from  the 
dies  after  forging.  If  the  part  geometry  does  not  permit 
ejection,  modifications  of  part  geometry  have  to  be  done  to 
permit  ejection.  For  example,  a  part  with  undercuts,  reverse 
tapers  etc.  cannot  be  ejected.  For  purposes  of  forging,  this 
is  removed  by  appropiate  modifications,  and  then  handled  by 
subsequent  machining. 

Guide  lines  for  the  choice  of  forging  direction  are  as 
follows : 

(i)  The  part  must  be  ejectable  from  the  dies  after  forging. 

(ii)  The  plan  area  of  the  part  in  the  forging  plane,  which 

is  del ined  as  any  plane  normal  to  the  forging  direction,  is 
the  largest. 

(iii)  1  he  maximum  part  height  in  the  forging  direction  be 
within  the  limits  suggested  by  the  press  manufacturer.  It  is 
recommended  to  be  less  than  0.8  times  the  press  stroke. 

4.3.2  Region  Definition 

To  illustrate  the  definition  of  regions,  consider  the 
simple  part  configurations  in  Figure  20.  The  axisymmetric 
hub-flange  as  well  as  the  rib-web  part  have  the  same 
cross-section.  The  terms  rib-web  and  hub-flange  are 
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identified  with  respect  to  a  forging  direction  and  the 
geometry.  These  terms  may  correspondingly  be  defined  by 
profile  changes  in  the  cross-section.  The  demarkation  is 
approximately  shown  by  the  broken  lines  in  Figure  20. 

For  a  part  that  is  powder  forged,  the  part 
simplification  is  termed  "region  subdivision".  The  terms 
rib-web,  hub-flange  lose  significance  and  are  called  as 
regions.  Hence,  for  a  cross-section,  the  region  subdivision 
criterion  is  based  on  significant  changes  in  the  profile. 

In  the  present  work,  the  region  interfaces,  viz,  the 
surfaces  that  are  common  to  a  pair  of  adjacent  regions,  are 
approximated  as  planes.  These  are  oriented  such  that  any 
inter-region  metal  flow,  that  is  postulated,  occurs 
approximately  in  a  direction  normal  to  the  interface. 

Ihe  plane  of  the  cross-section,  called  the  definition 
plane,  may  have  two  orientations  with  respect  to  the  forging 
plane.  Firstly,  it  may  be  normal  to  the  forging  plane. 

Figure  20a.  In  this  case,  the  interfaces  may  be 
unambiguously  determined  to  be  parallel  to  the  forging 
direction  at  locations  of  significant  profile  change. 
Secondly,  the  definition  plane  may  be  parallel  to  the  forging 
plane.  In  this  case,  Figure  20b,  two  possible  interface 
combinations  are  indicated.  Either  choice  is  permissible; 
but  the  choice  will  affect  the  preform  design,  as  will  be 
elaborated  in  subsequent  sections. 
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4.3.3  Generic  Metal  Flow  Patterns 

Various  combinations  of  metal  flow  across  interfaces, 
together  with  the  specified  forging  direction  may  be  used  to 
define  generic  metal  flow  patterns.  These  patterns  are 
equivalent  to  basic  deformation  modes.  They  are  illustrated 
in  Figure  21. 

Figures  21a  and  21b  represent  a  back-extrusion  and 
forward  extrusion  mode  respectively.  In  either  case,  the 
interface  metal  flow  is  into  the  region,  with  the  preform 
deformation  is  in  a  direction  opposite  to  or  in  the  same 
direction  as  the  ram  motion.  In  Figure  21c,  the  interface 
flow  is  reversed,  and  we  have  a  lateral  extrusion  mode.  Note 
that  in  Figures  21a-21c,  the  definition  plane  is  normal  to 
the  forging  plane. 

In  Figure  21d,  the  definition  plane  is  parallel  to  the 
forging  plane.  The  inter-region  metal  flow  shown  represents 
a  lateral  extrusion.  The  flow  direction  is  inconsequential, 
and  if  reversed  will  also  result  as  a  lateral  extrusion.  To 
contrast  with  the  lateral  extrusion  mode  of  Figure  21c,  the 
orientation  of  the  definition  plane  with  respect  to  the 
forging  plane  is  used.  In  Figure  21c,  the  definition  plane 
is  normal  to  the  forging  plane,  hence  the  mode  is  called  the 
normal  lateral  extrusion  mode.  In  Figure  21d,  the  definition 
plane  is  parallel  to  the  forging  plane,  hence  the  mode  is 
called  the  parallel  lateral  extrusion  mode. 
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Figure  21.  Generic  metai  flow  patterns  in  powder  forging. 
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In  Figure  21e,  there  is  effectively  no  interface  metal 
flow.  Deformation  is  constrained  to  occur  within  each 
region.  Such  a  mode  is  termed  the  parallel  upset  mode. 

Figure  21f,  analogously  represents  the  normal  upset  mode, 
where  in  this  case  the  preform  is  free  to  deform  in  a 
direction  normal  to  the  plane  of  the  paper.  This  is  possible 
only  for  a  non-axisymmetr ic  part.  For  axial  symmetry,  the 
configuration  will  result  in  a  repress  deformation. 

The  metal  flow  patterns  which  occur  in  various  regions 
depend  upon  profile  differences  between  the  preform  with 
respect  to  the  rigid  tooling.  In  regions  where  clearances 
exist  between  the  preform  and  tooling  profiles,  metal  flow 
must  be  present  no  that  the  clearances  are  closed  during 
forging.  Hence,  clearances  may  be  explicitly  defined  to 
influence  the  metal  flow  during  forging. 

In  extrusion  modes,  inter-region  metal  flov;  is  present 
across  interfaces.  The  extent  of  flow  is  dictated  by  the 
mass  deficiencies  that  have  to  be  accounted  for.  On  the 
other  hand,  in  upset  deformation  modes,  there  is  no 
substantial  metal  flow  across  interfaces.  Hence,  all 
deformation  is  confined  within  the  region.  The  mass  of 
material  in  the  preform  for  the  upset  region  will  correspond 
to  the  finished  mass  in  the  region. 
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4.3.4  Assessment  of  Fracture  Potential 

For  the  generic  deformation  modes,  a  qualitative 
assessment  of  fracture  likelihood  and  their  potential  sites 
may  be  done  to  determine  optimum  preform  shapes.  These 
optimum  shapes  may  then  be  quantified  with  process 
parameters,  via  experimental  modeling  techniques,  to  develop 
guide  lines.  This  approach  was  pursued  in  detail  for  back 
extrusion  of  hub-flange  and  cup  configurations.  (3,4) 

We  now  restrict  our  attention  specifically  to 
non-ax isymmetr ic  shapes  that  exhibit  the  deformation  modes, 
as  shown  in  Figure  21.  Non-axisymmetr ic  shapes  are  best 
formed  by  restricting  metal  flow  to  one  direction. 
Unrestricted  flow  results  in  early  fracture,  due  to  the  poor 
workability  of  sintered  powder  preforms.  Restricting  flow 
provides  for  some  compression  effects,  which  inhibit  fracture 
initiation. 

For  the  deformation  modes  in  Figures  21a-21e,  metal  flow 
is  permitted  in  the  definition  plane.  Consequently,  flow 
constraint  is  imposed  in  a  direction  normal  to  the  definition 
plane.  Conversely,  in  Figure  21f,  there  is  no  flow  possible 
in  the  definition  plane,  and  hence,  flow  is  permitted  to 
occur  in  a  direction  normal  to  the  definition  plane.  If  no 
flow  is  explicitly  permitted  in  any  direction  in  Figure  21f, 
a  repress  deformation  is  obtained  as  a  special  case. 


In  back  and  forward  extrusion  modes,  the  expanding  free 
surface  is  susceptible  to  fracture  due  to  friction  at  the 
side  walls.  For  no-draft  extrusion,  there  is  no  compression 
generated  by  the  side  walls  to  inhibit  the  possible  growth  of 
tensile  stresses. 

The  vicinity  of  the  fillet  undergoes  rapid  changes  in 
deformation  intensity  as  the  material  moves  from  under  the 
punch  into  the  rib,  where  it  does  not  undergo  any  futher 
plastic  deformation.  Die  contact  surface  fracture  is  a 
possibility  if  the  radius  is  small. 

For  the  cases  of  normal  and  parallel  lateral  extrusion, 
die  contact  surface  fracture  is  less  probable  due  to  the 
presence  of  compression  in  the  fillet.  Free  surface  fracture 
is  also  suppressed  at  the  expanding  surfaces  due  to  the 
presence  of  compression.  Intrusion  defect  formation  is 
unlikely  for  these  modes. 

For  the  upset  modes,  there  is  no  appreciable  flow  in  the 
vicinity  of  the  fillet  to  initiate  die  contact  surface 
fracture.  Free  surface  fracture  is  the  only  type  of  fracture 
that  is  likely. 


4.4  Experimental  Criteria  for  Fracture  Prevention 

Of  all  the  generic  modes  of  deformation  encountered, 
back  extrusion  is  likely  to  exhibit  greater  fracture 


problems.  The  influence  of  process  parameters  and  preform 
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options  were  studied  in  a  parallel  experimental 
(3233) 

program.  •  '  Relavent  information  is  summarized  in  the 

following  paragraphs. 

4.4.1  Experimental  Conditions 

The  experimental  process  conditions  for  sintered  powder 
forging  of  single  rib  and  opposed  rib  configurations  are 
given  in  Table  1.  The  material  was  aluminum  alloy  601AB  and 
601AC  of  initial  density  approximately  85%  theoretical. 

Copper  dispersed  in  oil  was  the  lubricant. 

The  deformation  was  done  cold.  The  preforms  were  0.75 
inches  square  and  flat.  Each  preform  was  split  at  a  central 
plane  and  grid  marks  were  etched.  The  two  halves  were 
assembled  and  forged,  thus  providing  a  deformed  grid  pattern. 

4.4.2  Summary  of  Results 

Free  surface  fracture  occurs  during  very  early  stages  of 
rib  formation  in  the  forging  of  both  single  rib  and  opposed 
rib  shapes.  The  occurrence  was  not  affected  to  any  great 
degree  by  the  size  of  the  fillet  radii  or  the  initial  preform 
aspect  ratio.  The  %  reduction  in  area,  i.  e.,  the  rib 
thickness,  was  the  only  significant  parameter  which  affects 
to  any  noticable  degree.  The  trend  is  towards  earlier 
fracture  with  increasing  %  reduction  in  area  (thinner  ribs). 
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Table  1.  Experimental  process  conditions  for  single  and  opposed 
rib  forging. 
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Material  :  601AB,  601AC  aluminum  alloy  powder. 

Preform  density  :  85%  theoretical. 

Lubricant  :  Copper  dispersed  in  oil. 

Temperature  :  Room  temperature. 

Preform  aspect  ratios  (H0/w0)  :  0.6,  0.5>  0.*4,  0.3,  0.2 
%  Reduction  in  area  (%R)  :  30%,  60%,  90% 

Fillet  radius  (r)  :  sharp,  1/64",  1/16",  1/8" 
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Die  contact  surface  fracture,  as  expected,  was  primarily 
controlled  by  the  punch  corner  radius.  Smaller  radii 
promoted  the  occurrence  of  fracture.  The  l/8th  inch  radius 
case  did  not  exhibit  any  die  contact  surface  fracture. 

The  preform  aspect  ratio  indicated  some  correlation  to 
the  occurrence  of  side  intrusion  defects.  In  general,  the 
occurrence  was  increased  with  decreasing  preform  aspect 
ratios  and  %  reduction  in  area.  For  single  rib  shapes, 
bottom  intrusion  occurs  for  low  preform  aspect  ratios  but 
higher  %  reduction  in  area.  For  opposed  rib  shapes,  these 
came  conditions  influence  internal  fracture.  The  scatter  band 
of  the  data  was  too  wide  to  derive  any  quantitative 
correlations . 

The  rather  premature  occurrence  of  free  surface  fracture 
led  to  a  closer  examination  of  the  growth  of  strain  at  the 
surface  during  deformation.  The  strains  were  measured  and 
are  plotted  in  Figure  22,  for  a  thin,  single  rib. (33) 

For  the  two  punch  corner  radii  indicated,  the  strain 
increases  very  rapidly  and  tapers  off  at  larger  punch 
displacements.  The  workability  limit  for  the  material  is 
shown,  and  it  is  evident  as  to  why  fracture  was  exhibited. 

The  rapid  growth  and  subsequent  tapering  off  suggests  a 
corner  radius  influence.  Specifically,  that  as  the  material  is 
rounding  the  corner,  the  free  surface  undergoes  excessive 
strain.  To  control  the  initial  strain  effects,  a  preform 
with  the  corner  built  in  is  inferred.  Such  a  preform  was 
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during  forging. 


subsequently  used,  and  the  results  are  shown.  The  importance 
of  preform  design  is  once  again  demonstrated.  Note  that  the 
ribs  are  draftless.  The  presence  of  a  draft  angle  increases 
the  influence  of  friction  and  the  strain  at  the  free  surface 
will  continously  increase  with  increasing  deformation. 

Densif ication  of  thin  ribs  was  good,  due  to  the  large 
amounts  of  lateral  flow  at  the  rib  base.  As  the  thickness 
increases,  greater  density  gradients  are  to  be  expected 
across  the  thickness  of  the  rib.  This  implies  that  thick 
ribs  have  to  undergo  some  amounts  of  repressing  to  fully 
density  it.  Consequently,  the  forging  of  thick  ribs,  if 
possible,  is  to  be  preferably  done  by  upsetting. 

4.5  Overview  of  CAD  of  Preforms 

The  computer-aided  design  of  preforms  to  powder  forge  a 
complex  part  consists  of  two  major  activities,  namely,  (i) 
Part  Geometry  Description  and  (ii)  Preform  Design.  The  two 
may  be  independent  of  each  other,  however,  for  a  coherent 
transition  from  one  to  another  requires  that  the  part 
description  be  compatible  and  complete  to  the  extent  required 
by  the  preform  design  phase.  A  general  part  description  is 
not  the  goal  of  this  research.  Any  part  description 
simplifications  that  result  due  to  limitations  in  the  preform 
design  stage  are  expolited.  Simplicity,  without  loss  of 
accuracy,  was  the  emphasis  for  the  design  of  the  part 
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description  program. 

Conceptually,  the  preform  design  phase  may  be  roughly 
classified  to  consist  of  the  following  subtasks: 

(1)  Sub-division  of  a  complex  part  into  regions. 

(2)  Identification  of  basic  deformation  modes. 

(3)  Determination  of  preform  shapes  for  each  mode. 

(4)  Assembly  of  individual  preform  shapes. 

The  sub-division  of  a  complex  part  into  regions  may  be 
accomplished  either  by  the  designer  explicitly  specifying  the 
regions  or  the  program  determining  the  region  sub-divisions 
automatically  according  to  some  set  of  rules.  At  this  stage, 
there  is  not  enough  experience  to  be  able  to  define 
completely  the  region  sub-division  rules  for  the  latcer 
process.  Consequently,  region  sub-division  is  delegated  to 
the  designer.  The  regions  are  explicitly  soecified,  during 
part  description,  according  to  guidelines  which  need  to  be 
adhered  to,  to  prevent  possible  design  failures  at  a  later 
stage . 

Assuming  that  the  complex  part  has  been  sub-divided  into 
regions,  two  routes  may  be  stipulated  towards  the 
determination  of  the  preform  shape.  The  first  route  is  the 
automatic  preform  design,  wherein,  generic  deformation  modes 
are  identified  according  to  the  region  configuration;  and  for 
each  mode  a  specific  preform  shape,  that  is  established  as 
optimum,  is  automatically  determined.  This  process  would 
require  no  user  intervention.  The  complexity  and  the 
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relative  inexperience  in  the  vast  possibilities  preclude  such 
an  ideal  scheme.  The  second  route,  is  the  interactive 
approach  approach,  where  the  preform  shape  is  arbitrarily 
specified  for  each  region  by  the  designer.  Depending  upon 
the  specification,  the  generic  modes  are  isolated  and  the 
speci 1 ieo  shape  is  evaluated  for  consistency,  densif ication, 
fracture  etc.,  by  the  program.  The  preform  profile  is  then 
interactively  manipulated  by  the  designer  till  a  satisfactory 
preform  is  determined.  Some  experience  is  required  to  judge 
whether  the  preform  is  satisfactory  from  a  practical 
standpoint.  However,  in  this  approach,  additional  decision 
making  functions  can  be  built  in  as  further  knowledge  and 
experience  is  obtained  in  preform  design.  The  accuracy  of 

the  program  will  have  to  evolve  through  constant  additions  in 
the  future. 

Figure  23  illustrates  the  overall  structure  of  the 
present  CAD  program.  The  part  description  phase  is  a  seperate 
activity,  primarily  because  of  the  probability  of  error  being 
higher.  The  preform  design  phase  is  in  an  interactive 
graphics  environment  to  allow  for  a  rapid  evaluation  of  trial 
preforms.  The  details  are  now  considered  in  the  appropriate 


sequence . 
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Figure  23.  Overall  structure  of  computer-aided  preform  design  program. 
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4.6  Part  Geometry  Description 

A  part  geometry  description  program  was  developed  to 
describe  a  3-dimensional  part  in  a  form  that  is  suitable  for 
subsequent  processing  by  the  preform  design  phase.  The 
program  emphasizes  simplicity  and,  hence,  limitations  have 
been  imposed  to  accomplish  the  description.  The  choice  of 
APT  for  part  description  was  rejected.  Firstly,  because  it 
is  designed  for  use  in  machining,  and  an  APT  description 
would  require  considerable  reprocessing  to  convert  to  a  form 
that  is  suitable  for  forging  preform  design.  Secondly,  it 
needs  considerable  user  training.  Thirdly,  the  size  of  the 
APT  program  is  much  too  large  as  it  contains  features 
completely  unrelated  to  forging. 

The  program  that  was  developed  is  called  "PADEL"  (PArt 
DEscription  Language).  It  is  tailored  to  specifically  meet 
the  requirements  of  powder  forging.  In  addition  to  being 
able  to  describe  the  part,  region  sub-divisions  may  be 
readily  accomplished,  forging  directions  specified,  volumes 
calculated  etc.,  which  are  features  peculiar  to  powder 
forging.  It  is  not  generalized  to  describe  all  possible 
components.  Nevertheless,  it  is  capable  of  describing  a 
large  class  of  parts  that  are  currently  being  powder  forged. 


4.G.1  Description  Space 
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The  part  description  is  done  in  the  first  octant  of  a 
3-dimensional  Carthesian  space,  called  the  description  space, 
gure  24.  The  coordinate  planes  are  referred  to  as  XYPLN, 
YZPLN  and  ZXPLN.  The  forging  plane  is  implicitly  fixed  as 

XYPLN,  hence,  the  forging  direction  is  the  negative 
Z-direction. 

l>rior  to  description,  the  part  is  oriented  so  that  it  is 
wholly  i„  the  description  space  and  has  the  desired  forging 
direction.  If  the  part  has  any  plane  of  symmetry,  then  the 
appropriate  plane  is  made  to  coincide  with  a  coordinate 
Plane.  The  coordinate  plane  is  then  explicitly  identified  as 
being  a  plane  of  symmetry.  In  this  event,  only  that  portion 
of  the  part  in  the  description  space  need  be  defined. 

i.6.2  Zone  Description 


A  ^no  is  oe fined  as  c,  volume  th£|t  is  gonerated  by 

translating  or  rotating  about  an  axis  a  representative 

ss  section.  The  translation  or  rotation  is  bounded  by  two 
planes,  called  the  primary  and  secondary  definition  planes. 
Hence,  in  effect,  a  zone  is  describable  by  defining  the 
or imaLy  and  secondary  definition  planes,  along  with  the 
cross-section  profile  in  the  primary  definition  plane. 

Any  part  that  is  describable  by  PADETL  is  assumed  to  be 
decomposable  into  zones.  The  decomposition  of  a  part  in 
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Figure  2k.  Definition  of  the  part  description  space. 
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Figure  25. 


Primary  definition  plane  orientations  possible  for  the 
definition  of  zones. 
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description  space  is  accomplished  by  defining  "cut  planes" 
that  intersect  the  part,  thereby  isolating  the  zones.  The 
cut  planes  are  always  normal  to  the  forging  plane,  and 
furthermore,  do  not  permit  any  metal  flow  from  one  zone  to 
another  through  itself.  The  pair  of  definition  planes  that 
define  any  zone  may  be  either  parallel  or  normal  to  the 
forging  plane,  Figure  25. 

By  definition  of  a  zone,  a  restriction  is  placed  on  the 
geometric  complexity  of  the  part  that  can  be  described. 
Specifically,  variation  in  the  cross-section  profile  during 
translation  or  rotation  cannot  be  represented.  The  advantage 
of  this  restriction  is  that  a  complex  3-dimensional  part  may 
be  considered  as  a  series  of  cross-sections.  Preform  design 
is  considerably  simplified  when  done  on  a  cross-section 
level.  However,  small  features  in  the  part,  which  cannot  be 
decomposed  into  zones,  can  be  indirectly  input  as  an  excess 
volume  to  be  added  on  to  the  appropriate  zone. 


4.6.3  Cross-Section  Description 


Any  zone  is  completely  defined  by  its  definition  planes 
and  the  cross-section  profile  in  the  primary  definition 
plane.  The  cross-section  profile  is  described  by  two 
contours,  whose  ends  are  connected  by  straight  lines  to  close 
the  profile.  Such  a  description  serves  to  distinguish  the 
surfaces  of  the  zone  that  contact  the  punches,  die  side  walls 
and  adjacent  zones. 
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In  the  case  where  the  primary  definition  plane  is  normal 
to  the  forging  plane,  the  contours  that  are  described 
represent  intersections  of  the  primary  definition  plane  with 
the  upper  and  lower  punch  surfaces.  The  ends  of  the  contour, 
which  are  connected  by  straight  lines,  in  turn,  represent  the 
intersection  of  the  definition  plane  with  either  the  die  wall 
surfaces  or  cut  planes. 

In  the  case  where  the  primary  definition  plane  is 
parallel  to  the  forging  plane,  the  contours  represent 
intersection  of  the  primary  plane  with  the  die  wall  surfaces. 
The  cross-sections  on  the  definition  plane,  in  themselves, 
are  surfaces  that  are  coincident  with  the  upper  and  lower 
punches  respectively. 

Each  of  the  two  contours  are  described  using  statements 
borrowed  from  APT.  The  contour  is  assumed  to  be  a  contiguous 
sequence  of  line  segments  and  circular  arcs.  Each  of  these 
segments  are  first  explicitly  defined  by  the  respective  APT 
geometric  statements  for  lines  and  circles  and  then  the 
sequence  is  specified  by  other  statements. 

4.6.4  Region  Description 

Once  the  two  contours  of  a  cross-section  are  defined, 
the  region  sub-division  is  accomplished  according  to  changes 
in  the  contour.  The  sub-division  is  done  by  specifying  the 
interface  locations.  The  interfaces  are  strictly  plane 
surfaces,  which  upon  intersection  with  the  primary  definition 
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plane,  result  in  a  line  that  intersects  the  two  contours  only 
once.  The  interface  lines  are  defined  so  that  the  lateral 
metal  flow  is  approximately  normal  to  the  interface  line. 

'.6.5  Structure  of  PADEL 

The  part  description  program,  PADEL,  is  implemented  as  a 
two  stage  processor.  The  input  to  the  processor  is  a  source 
program  consisting  of  legal  statements  (Appendix  C) ,  that 
describe  the  part. 

Stage  I  processing  is  essentially  a  syntax  check  on  the 
validity  of  the  source  program  statements.  Each  statement 
has  a  rigid  syntax  that  must  be  adhered  to.  Any  syntax 
errors  detected  during  processing  will  prevent  entry  into 
stage  II. 

Stage  II  processing  is  the  arithmetic  processing  stage. 
The  part  is  internally  reconstructed  by  the  program  as  an 
assembly  of  regions.  The  surfaces  of  each  region  are 
identified  with  punch  surfaces,  die  wall  surfaces,  cut  planes 
etc.  Any  errors  detected  at  this  stage  will  cause  PADEL  to 
abort.  On  successful  processing,  the  part  is  stored  on  a 
zone  by  zone  basis  in  a  disk  file  for  use  by 
design  phase. 


the  preform 
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4.7  Preform  Design  Methodology 

Preform  design  is  accomplished  on  a  region  by  region 
basis.  The  flow  of  metal  in  the  regions  is  determined  by  the 
relative  shape  difference  between  the  tooling  and  the 
preform.  Tooling  profiles  are  established  in  the  part 
description  stage.  Hence,  once  an  approximate  preform  shape 
is  provided,  metal  flow  may  be  inferred  and  evaluated. 

The  description  by  PADEL  reduces  a  complex  component  to 
a  set  <'f  zones.  The  cross-section  of  each  zone  consists  of 
two  contours,  called  the  TOP  and  BOTtom  contour  respectively, 
which  are  connected  at  its  end  points  by  straight  lines. 

Each  zone  is  further  decomposed  into  regions,  the  interfaces 
being  represented  in  the  cross-section  plane  by  interface 
lines  that  intersect  the  contour  at  points.  These 
intersection  points  between  the  region  interface  lines  ana 
the  contours  are  called  "nodes".  Note  that  the  end  points  of 
the  two  contours  also  become  nodes.  By  virtue  of  these 
definitions,  the  number  of  nodes  on  either  contour  is  one 
more  than  the  number  of  regions  defined  in  the  zone;  the 
number  of  nodes  on  both  contour  being  equal.  Figure  26 
illustrates  these  definitions. 

The  specification  of  a  preform  shape  is  a  user  function, 
done  in  an  interactive  graphics  environment.  Recall  that  the 
contours  of  a  cross-section  represent  the  intersection  of  the 
primary  definition  plane  with  the  tooling  surfaces  at  the  end 
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FILL  /TOP, 21, FROM, 2, 8.0 
FILL  /  TOP, 34, FROM, 3, 5.0 
FILL  /  BOT, 34, FROM, 3, 5.0 


Figure  26.  Definition  of  nodes  on  the  'TOP'  and  'BOT'  contour  of  a 

zone  cross-section.  The  result  of  a  preform  specification 
by  the  user,  with  respect  to  the  contours,  is  shown. 


\ 
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of  the  forging  stroke.  The  preform  shape  input  is  done  with 
respect  to  the  contours,  and  in  effect,  represents  the 
contact  configuration  of  the  tooling  with  the  preform  at  the 
start  of  the  forging.  Preform  design  then  entails  a  proper 
seperation  of  the  user  input  preform  shape  profiles,  so  that 
proper  weight  control  and  distribution  is  obtained,  and  the 
evaluation  for  non-aggravating  metal  flow  and  full 
densi f ication . 

4.7.1  Preform  Shape  Description 

The  contours  of  a  zone  are  defined  in  a  local  x-y  plane. 
The  nodes  are  ordered  left  to  right  by  increasing  values  of 
its  x-coord inates .  Between  a  pair  of  nodes  is  the  tooling 
profile  for  the  region(s).  Figure  26. 

Preform  shape  description  may  be  accomplished  in  two 
directions.  The  first  is  the  "filling  of  cavities"  of  the 
top  and  bottom  contours.  The  second  is  the  specification  of 
a  clearance  with  respect  to  the  end  lines  of  the  contours. 

The  terms  'top'  and  'bottom'  are  relative,  and  are  defined 
according  to  the  physical  format  of  the  display  on  the  CRT. 

It  is  recommended  that  they  correspond  to  the  top  and  bottom 
punches.  The  present  version  of  the  preform  shape  evaluation 
program  assumes  this  correspondence. 
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Consider  an  arbitrary  cross-section,  each  of  whose 
contours  have  M  nodes  ordered  in  increasing  x-coordinates . 

(i)  The  level  of  any  node  is  defined  as  the  value  of  the 
y-coordinate  of  the  node,  with  respect  to  the  local  reference 
frame . 

(ii)  The  minimum  base  level  of  a  contour  is  defined 
according  to 


(4-1) 


(iii)  Any  profile  between  two  nodes  on  either  contour  is 
referenced  by  the  ordered  node  pair  IJ;  where  strictly,  I  * 

J  . 

(iv)  Any  interface  line  that  connects  nodes  on  the  top  and 
bottom  contours  is  referenced  by  the  node  pair  IJ,  where 
strictly,  I  =  J. 

(v)  A  cavity  is  said  to  exist  between  the  nodes  I  and  J  of  a 
contour  if  the  maximum  (minimum)  of  the  profile,  yp ,  is 


— TOP  Contour 


(4-2) 


y  <  v,  — BOT  Contour 
Jp  jrB 


P 


where  v  is  the  current  base  level  for  the  profile.  The 


cavity  depth,  Ii ,  is  measured  from  the  base  level  by 


(4-3) 


93 


Kith  these  definitions  the  preform  shape  input  scheme 
may  be  introduced.  The  shape  input  is  based  on  appropriate 
manipulations  of  the  base  level  of  each  contour,  from  the 
minimum  level,  by  means  of  explicit  commands.  By 
progressively  updating  the  base  level,  various  preform  shapes 
may  be  input.  Three  commands  are  available,  namely,  the  FILL/ 
command,  the  FLOAT/  command  and  the  NOFLOW/  command.  Their 
description  and  function  follow. 

Filling  of  cavities  is  accomplished  by  a  command  of  the 

form 


FILL/ [ TOP, BOT] ,IJ,FROM, l I , J ] ,k 
where  [TOP, BOT]  represent  an  explicit  identification  for 

the  top  and  bottom  contour  respectively. 

IJ  references  the  node  pair,  I=J. 

FROM, (I, J]  explicitly  changes  the  base  level 
to  the  referenced  node  I  or  J. 
k  is  the  degree  of  fill,  measured 

from  the  base  level,  that  is  desired 
and  0<k<10.  A  unit  degree  of  fill 
is  arbitrarily  defined  as  H/10. 

Specification  of  k=0  will  result  in  the  preform  shape 
between  the  nodes  I  and  J  to  be  a  horizontal  line  through  the 
reference  node.  For  k=10,  the  preform  shape  matches  the 
tooling  profile,  i.e.,  no  clearance  is  provided. 

Intermediate  values  of  k  will  result  in  a  clearance. 
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Once  a  profile  is  'filled'  by  the  command,  it  can  only 
be  'refilled'  or  'unfilled'  by  restarting  the  shape  input. 
However,  other  profiles  of  a  contour  may  be  addressed  and 
manipulated.  The  current  base  level  is  used  to  automatically 
define  the  preform  shape  for  those  profiles  not  addressed  by 
an  explicit  command.  Figure  26  illustrates  the  result  of  a 
FILL/  specification.  Note  that  the  FILL/  command  is 
exclusively  used  to  describe  the  preform  shape  with  respect 
to  the  top  and  bottom  contours  only. 

Floating  of  the  preform  from  the  end  lines  of  the 
contours  is  accomplished  by  a  command  of  the  form 

FLOAT/I I, k 

where  II  represents  the  node  pair  corresponding 
to  the  end  line. 

k  is  the  degree  of  float,  measured 

from  the  end  line,  that  is  desired  and, 

0<k<10.  A  unit  degree  of  float  is 
defined  arbitrarily  as  1/1 0-th  the 
distance  from  the  end  line  to  the 
interface  line. 

Specification  of  k=0,  will  result  in  the  preform  being 
flush  with  the  end  line.  For  k=10,  the  preform  is  flush  with 
the  interface  line,  and  hence,  is  completely  absorbed  in  the 
adjacent  region. 
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To  provide  for  flow  in  a  direction  normal  to  the 
definition  plane,  the  command 

NOFLOW/ I I 

where  II  represents  an  interface  line,  is  used.  The  preform 
design  program  always  attempts  to  provide  flow  in  one 
direction.  Hence  by  properly  constraining  flow  in  the 
definition  plane  and  a  judicious  use  of  NOFLOW/,  flow  in  the 
normal  direction  may  be  achieved. 

The  use  of  the  above  commands  is  illustrated  in  a  later 
section  where  the  preform  shape  is  specified  for  a  weapon 
component . 

4.7.2  Preform  Shape  Evaluation 

The  evaluation  of  the  shape  input  by  the  user  is  done  in 
two  stages.  First,  the  stroke  is  computed  by  appropriate 
mass  balances  and  second,  the  overall  metal  flow  that  occurs 
in  the  regions  is  checked  for  consistency  with  predetermined 
allowable  conditions.  Inconsistent  metal  flow  will  lead  to 
undensified  regions.  These  situations  are  rectified  by  the 
user  by  appropriate  manipulations  on  the  input  shape. 

The  bottom  punch  is  assumed  to  be  stationary.  The 
initial  preform  density  is  assumed  to  be  constant  for  all 
regions  in  the  part.  The  program  attempts  to  determine  a 
single  stroke  for  the  upper  punch  that  is  required  to  size 
the  preform.  In  the  case  of  hot  repressing  of  multi-level 
parts  with  preforms  of  uniform  initial  density,  a  single 
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stroke  will  result  in  undensified  regions.  In  such  cases, 
split  punches  are  used  where  each  punch  has  a  stroke  that 
provides  for  full  densif ication  of  a  region.^)  a  single 
stroke  option  for  repressing  is  applicable  only  in  the  case 
where  the  preform  has  a  non-uniform  initial  density  in 
various  regions.  In  contrast,  presence  of  metal  flow  durinq 
upset  forging  of  preforms  allows  for  a  single  stroke  anj 

uniform  initial  preform  density.  This  necessitates  that 
I- roper  clearances  for  metal  flow  be  provided  so  that  a  proper 

preform  weight  distribution  is  achieved  to  form  a  fully 
densifieu  port. 

Any  region  of  the  part  consists  of  6  surfaces,  as  shown 
in  Figure  27.  Each  region  exhibits  a  characteristic  metal 
flow  during  forging,  according  to  the  relative  shape 
difference  (and  hence  clearance)  between  the  preform  and  the 
tooling.  Table  2  illustrates  the  various  metal  flow 
possibilities  that  occur  in  a  region  during  forging.  These 
flow  possibilities  are  for  regions  derived  from  a  parent 
zone,  whose  primary  definition  plane  is  normal  to  the  forging 
plane.  Consequently,  region  interfaces  are  coincident  with 
the  Left  and  Right  surfaces.  The  presence  of  metal  flow 
across  interfaces  is  used  to  classify  a  region  as  a  member  of 
an  extrusion  deformation  mode  (or  extrusion  group) .  No  metal 
flow  across  interfaces  classifies  a  region  as  a  member  of  an 
upset  deformation  mode  (or  upset  group).  Note  that  similar 
metal  flow  possibilities  exist  for  regions  derived  from  a 
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Table  2.  Metal  flow  directions  in  regions  during  forging. 


DUAL  EXTRUSION 
Interface  flow  into  region 


BACK  EXTRUSION 
Interface  flow  into  region 


FORWARD  EXTRUSION 


Interface  flow  into  region 
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LATERAL  (normal)  EXTRUSION 
Interface  flow  into  region 
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parent  zone  whose  primary  definition  plane  is  parallel  to  the 
forging  plane.  These  possibilities  are  not  incorporated  in 
the  present  shape  evaluation  program.  The  extension  is  a 
future  task,  and  may  be  readily  appended  by  following 
procedures  analogous  to  those  currently  implemented. 

The  mass  balance  strategy  for  a  region  is  dependent  upon 
whether  the  region  is  a  member  of  an  extrusion  group  or  an 
upset  group.  For  an  arbitrary  region,  the  mass  of  material 
in  the  region  at  the  finished  forged  state,  M  ,  is  given  by 

Mf  =  pf(Af*Df  +Pf)  (4-4) 

where  p  =  theoretical  density  of  the  material. 

A  =  cross-sectional  area  of  the  region. 

D  =  depth  of  region  in  the  normal  direction, 
r  =  additional  volume  for  the  region,  due  to  a  feature 
that  could  not  be  input  via  PADEL. 

The  mass  of  preform  material  required  in  the  region  need 
not  be  equal  to  M  .  This  is  true  when  the  region  belongs  to 
an  extrusion  group  and  exchanges  mass  across  interfaces. 
However,  if  it  belongs  to  an  upset  group,  there  is  no  mass 
exchange,  in  which  case  the  mass  of  preform  material  is  equal 
to  M  .  Recall  that  the  preform  shape  is  input  with  respect 
to  the  tooling  profiles  at  the  end  of  the  stroke.  Hence,  the 
cross-section  area  enclosed  by  the  preform  input  shape  for  a 
region,  A  ,  is  less  than  or  equal  to  Af .  This  depends  upon 
the  shape  specified  by  the  user.  Figure  28  illustrates  three 
possible  cases  of  shape  (clearance)  specification. 
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Figure  28.  Schematic  representat ion  of  stroke  computations  for  regions 
with  clearances  in  each  of  the  three  directions. 
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For  a  stroke, A,  the  mass  of  the  preform  material  in  a 
region  is  given  by 


M 


p^W  '  pp(A-Vdp 


(4-5) 


where  wp  is  the  width  as  shown  in  Figure  23.  The  subscript 

’p’  corresponds  to  quantities  defined  for  the  preform. 

Hence,  for  the  region  in  an  upset  group,  M  =Mf.  For  the 

region  in  an  extrusion  group,  M  =  M  .  Exchange  of  mass  is 

P  i 

postulated  to  occur  out  of  a  region  through  interfaces  if  M 

P 

>Mf  .  Conversely,  if  Mp  <Mf,  mass  exchange  is  into  the  region. 

A  region  that  is  a  member  of  an  extrusion  group  must  be 

paired  with  another  region  of  the  same  group,  across  the 

interface  that  permits  metal  flow.  The  direction  of  metal 

flow  is  from  a  region  of  excess  mass  to  a  region  of  deficit 

mass.  The  total  mass  of  preform  material  in  an  extrusion 

group,  must  equal  that  required  by  the  group  in  the  finished 

part.  This  mass  balance  is  given  by 

L  L  L 

p  E  (AD  +  T  )  =  p  {  E  (AD).  +  A  E  (w  D  ).  }  (4-6) 

£  i-K  f  f  f  1  P  i-K  P  P  1  f=K  P  P  1 

By  definition  of  the  group,  (D  ).  =(D  )  ,  and  the  stroke,  a, 

pi  f  i 

for  the  group  is 


=  {  E(v  )  .  /  D  -  E  (A  )  .  }  /  E(w  )  . 

f  i  pi  pi 


(4-7) 


where  p=p  /  p„  is  the  fractional  density  of  the  preform. 

p  f 

d=D^=D^  is  the  depth  of  the  extrusion  group. 

<Vf  <AfDf*r£>i 
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Once  the  stroke  is  computed,  tne  metal  flow  directions 
across  interfaces  may  be  readily  determined  by  computing  the 
excess  or  deficit  mass  for  each  region  from  (4-4)  and  (4-5). 
For  a  valid  preform,  the  direction  of  metal  flow  must  be 
consistent  with  the  allowable  flow  configurations  listed  in 
Table  2.  An  inconsistent  flow  pattern  will  imply  the 
possibility  of  under  densif icat ion ,  this  situation  is 

rectified  by  appropriate  modification  of  the  input  shape  by 
the  user. 

A  region  that  belongs  to  an  upset  group  has  no  mass 
exchange  with  an  adjacent  region.  The  required  mass  is  fully 
contained  in  the  preform.  The  stroke, A  ,  is  directly  related 

to  the  clearance.  The  mass  balance  for  the  region  is  given 
by 

pf(AfDf  +  rf)  =  pp{ApDp  +  (A.wp)Dp  }  (4-3) 

Hence,  the  stroke, a  ,  is 


A  =  {v  /  pd  -  A  }  /  w 

f  P  p  p 


(4-9) 


clearances  can  occur  in  two  directions,  as  shown  in 

Figure  28.  Firstly,  it  can  be  in  the  lateral  direction  and 

secondly,  it  can  be  in  a  direction  normal  to  the 

cross-section  plane.  In  the  first  case,  we  have  D  =D  =D,  and 

P  f 

the  clearance  'c*  is  given  by  c  =  (w  -  w  ).  The  stroke  is 

t  p 

then  given  by 


A 


(4-9a) 


102* 


=  (  V,/  pD  -  A  )  /  (wtf  -c) 
f  p  f 


where  A  is 
P 

In  the 

A  =A  .  The 
P  f 

stroke  is 


an  imDlicit  function  of  c. 

second  case,  we  have  =wf  =w»  and  consequently 

clearance  'c1  is  defined  by  c=(D  -D  ).  The 

f  P 


A  =  {v./  p(D.-c)  -  A  }  /  w  (4-9b) 

t  r  P  P 

In  both  equations  (4-9a)  and  (4-9b) ,  the  stroke  is 
dependent  upon  the  clearance.  We  can  define  a  stroke  range 
that  varies  between  a  minimum  value  corresponding  to  a  jero 
clearance  and  a  maximum  value  corresponding  to  a  maximum 
allowable  clearance.  The  maximum  allowable  clearance  is 
strictly  to  be  determined  from  workability  considerations. 
Presently  it  is  set  to  be  half  the  total  length  in  the 
direction  of  flow. 

Consider  a  zone,  which  is  an  assembly  of  regions  that 
belong  to  the  upset  group.  Each  region  of  the  zone  will 
yield  a  stroke  range  if  no  clearances  are  specified 
explicitly.  On  the  other  hand,  if  clearances  are  explicitly 
specified,  then  each  region  will  yield  a  unique  stroke 
according  to  (4-9a)  or  ( 4  —  9 b ) .  These  strokes  need  not  be  the 
same,  and  we  have  a  situation  that  is  found  in  repressing. 

If  instead  of  specifying  explicit  clearances,  suppose  only 
the  desired  clearance  directions  are  spec-ified.  Then  for 
each  region,  a  maximum  stroke  and  a  minimum  stroke  may  be 
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obtained.  If  the  stroke  for  the  zone  is  chosen  as  the 
largest  of  the  minimum  strokes,  then  only  the  region  whose 
minimum  stroke  corresponds  with  the  chosen  stroke  undergoes 
repress.  The  other  regions  have  clearances,  which  may  be 
determined  by  solving  (4-9a)  and  (4-9b) .  It  is  for  this 
reason  that  a  specification  of  k=-l  in  the  FLOAT/  command  is 
available.  In  such  cases,  the  program  computes  the 
clearances  for  the  regions. 

If  the  zone  has  an  extrusion  group,  the  stroke  is 
explicitly  computed  by  (4-7).  This  stroke  then  forms  the 
basis  for  calculating  the  clearances  for  all  the  upset  group 
regions. 

In  the  case  of  a  multi-zone  part  with  only  upset  group 
regions,  the  optimum  stroke  is  chosen  from  the  optimum  for 
each  zone.  If  two  zones  have  extrusion  groups  with  different 
strokes,  then  user  intervention  is  necessary  in  choosing  the 
stroke  and  the  appropriate  shape  modification  required  to 
obtain  a  single  stroke  value. 

During  preform  shape  input,  the  program  monitors  each 
profile  description  as  it  is  input  by  the  user  and  attempts 
to  define  and  classify  the  regions.  If  an  extrusion  group 
necomes  defined,  an  attempt  to  change  to  an  upset  group  by 
the  use  of  the  NOFLOW/  command  will  be  foiled.  Being  an 
interactive  program,  various  preform  shapes  may  be  readily 
input  and  evaluated.  Upon  successful  design,  the  preform  is 
stored  on  a  zone  by  zone  basis  for  subsequent  compaction. 
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4.8  Application  to  a  Weapon  Component 

The  computer-aided  design  approach  is  applied  to  a 
non-ax isymmetr ic  weapon  component.  This  serves  to  determine 
the  feasibility  of  the  approach  and  provide  the  framework  for 
improvement  of  the  interaction  between  the  designer  and  the 
program. 

The  part  drawing  for  the  component  is  shown  in  Figure 
29.  The  first  task  is  the  description  of  the  component, 
along  with  identification  of  a  forging  direction  and  a 
schematic  sub-division  into  zones  and  regions.  The  second 
task  is  the  design  of  a  preform  using  the  CAD  approach.  The 
third  and  final  task  is  the  trial  forging  of  the  preform(s) 
to  evaluate  for  occurrence  of  defects  and  undensified 
regions. 

4.8.1  Geometry  Description 

The  weapon  component  has  a  plane  of  symmetry.  The  part 
is  shown  in  isometric,  Figure  30,  properly  oriented  in  the 
description  space.  Note  that  the  plane  of  symmetry  is  made 
coincident  with  a  coordinate  plane,  hence  only  half  the  part 
need  be  described. 

The  orientation  of  the  part,  as  shown  in  Figure  30, 
implies  that  the  forging  direction  has  been  decided.  For 
this  component,  this  is  the  only  orientation  that  will  permit 
ejection  of  the  part  after  forging.  The  presence  of  a 
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Figure  29. 


Part  drawing  for  a  non-symmetr ic  weapon  component. 
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Half  isometric  view  of  the  part  oriented  in  description 
space.  Note  the  plane  of  symmetry,  and  hence  only  half 
the  part  is  described. 
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protrusion,  called  lug,  at  the  bottom  left  precludes  forging 
in  the  X-direction.  However  if,  for  the  current  orientation, 
a  suitable  preform  cannot  be  obtained  then  the  forging 
direction  may  be  changed  with  subsequent  machining  of  the 
lugs.  In  the  present  orientation,  only  the  hole  in  the 
center,  Figure  29,  needs  to  be  drilled. 

Figure  31  illustrates  the  decomposition  of  the  part  into 
zones  by  slicing  with  cut  planes,  CTPLl  and  CTPL2 .  The  cut 
planes  are  normal  to  the  forging  plane.  The  part  is  shown  in 
an  exploded  view,  with  the  zones  isolated,  in  Figure  32.  The 
definition  planes  and  the  local  coordinate  planes  are  easily 
identified. 

At  this  point,  it  is  evident  that  the  lug  cannot  be 
described  as  a  zone  since  it  does  not  have  a  characteristic 
cross-section.  Therefore,  the  part  is  modified  for 
description  by  removing  the  lug.  The  removed  volume  of  metal 
is  then  input  seperately,  to  zone  1. 

The  contours  of  each  zone,  Figure  33,  are  now  described. 
Cach  contour  is  assumed  to  consist  of  a  contiguous  sequence 
of  line  segments  and  circular  arcs.  These  are  described 
using  statements  to  define  the  lines  and  circles  and  then  the 
connectivity  is  established.  The  statements  are  similar  in 
format  to  those  of  APT,  and  are  listed  in  Appendix  C. 

The  PADEL  source  program  that  describes  the  weapon 
component  is  given  in  Appendix  D.  Figures  31,  32  and  33  used 
in  conjunction  with  the  source  program  will  clarify  the 
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Fi gure  31 . 


Decomposition  of  the  part  into  zones  using  'cut  planes'. 


FDliGrj.fr,  I1RECTIDN 


111 


CM 

CO 

<D 


=3 

cn 


Exploded  view  of  the  part  showing  the  zones  and  the  definition  planes. 
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Figure  33.  Geometric  elements  for  the  contour  description  of  the  zones 
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uescription  scheme.  Note  that  the  lines  in  the  source 
program  that  begin  with  a  semi-colon  are  comment  lines  which 
are  ignored  by  PADEL.  The  results  of  processing  are  stored  in 
a  disk  file  for  input  to  the  preform  design  phase.  Graphics 
verification  of  the  input  is  available. 

4.8.2  Preform  Design 

The  output  from  PADEL  is  the  input  for  the  design  phase. 
Each  zone  of  the  part  is  explicitly  called  by  the  user, 
preform  shapes  are  input  and  evaluated.  Final  assembly  of 
the  individual  preforms  to  provide  the  overall  preform  is 
then  automatically  accomplished. 

When  a  zone  is  explicitly  called  for  preform  shape 
input,  the  top  and  bottom  contours  are  displyed  on  the  CRT 
along  with  the  nodes  and  the  schematic  region  sub-sivision  as 
shown  if  Figure  34.  The  preform  shape  input  commands,  to  be 
entered  by  the  user,  are  also  indicated. 

Consider  zone  1.  This  zone  is  divideo  into  3  regions, 
which  are  ordered  left  to  right.  The  minimum  base  levels  for 
the  two  contours  are  computed.  For  the  top  contour,  the 
minimum  base  level  corresponds  to  the  level  of  nodes  2  or  3. 
For  the  bottom  contour,  the  minimum  base  level  corresponds  to 
the  level  of  node  1.  Filling  of  the  shape  is  done  with 
respect  to  these  defined  base  levels.  If  no  commands  are 
specif  ieci,  the  base  levels  are  used  to  compute  the  preform 
shape.  Specifically,  for  the  top  contour,  the  preform  shape 
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Figure  3^.  Shape  input  commands  for  the  design  of  the  preform  for 
the  part  under  consideration. 
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in  region  1  a  horizontal  line  through  node  2;  follows  the 
contour  from  node  2  to  node  3  in  region  2  and  is  a  horizontal 
line  at  level  of  node  3  in  region  3. For  the  bottom  contour, 
the  shape  is  a  horizontal  line  through  node  1  for  all  the 
regions  as  the  base  level  does  not  intersect  the  contour.  The 
region  interfaces  22  and  33  are  assumed  to  permit  flow  unless 
explicitly  stated  otherwise  with  a  NOFLOW/  command. 

The  first  FILL/  command  of  Figure  34  addresses  the  TOP 
profile  of  region  1.  Note  that  the  fill  reference  level  node 
is  now  node  1,  hence  the  base  level  is  changed  to  the  level 
of  node  1.  If  the  shape  is  to  match  the  contour  in  region  3, 
there  is  no  necessity  to  specify  another  riT, T/  co .  This 
is  true  only  if  the  contour  in  region  3  is  lov.er  than  the 
established  base  level.  The  FILL/  command  for  the  bottom 
contour  has  the  same  effect.  The  result  of  this 
specification  is  shown  shaded  in  Figure  34.  The  result  of  a 
different  specification  for  the  same  zone  is  illustrated  in 
Figure  2G. 

Suppose  that  the  NOFLOW/  command  is  not  given.  Then, 
the  regions  of  zone  1  constitute  an  extrusion  group.  The 
stroke  is  calculated  from  a  mass  balance.  The  calculated  mass 
is  compared  with  the  mass  required  in  each  region  in  the 
finished  part.  For  an  initial  preform  density  of  80% 
theoretical,  the  calculated  stroke  is  0.15  inches.  The 
computed  masses  for  the  regions  are: 
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Region  1  Region  2  Region  3 


Required  mass 


0.2202 


0.1568 


0.1947 


Preform  mass 


0.2172 


0.2028 


0.1517 


Di f f erence 


-0.0030 


0.0460 


-0.0430 


The  interface  metal  flow  is  from  region  2  to  both 
regions  1  and  3.  However,  most  of  the  flew  is  from  region  2 
to  region  3,  with  very  little  flow  from  region  2  to  region  1. 
Region  1,  then,  is  approximately  under  repress  with  about  20% 
height  reduction.  For  this  stroke  it  is  improbable  that 
region  3  will  be  fully  densified. 

Specifying  the  NOFLOW/  commands  to  the  above  preform 
changes  the  mode  of  deformation.  As  a  result,  metal  flow  in 
each  region  is  contained  and,  the  direction  of  flow  is  normal 
to  the  plane  of  the  paper.  bach  region  then  becomes  a  member 
of  the  upset  group.  For  an  initial  preform  density  of  80% 
theoretical,  the  stroke  is  computed  to  be  0.48  inches  and 
this  stroke  corresponds  to  the  largest  of  the  minimum  strokes 
for  the  regions.  Furthermore,  it  corresponds  to  the  minimum 
stroke  of  region  3.  Hence,  now  region  3  undergoes  a  repress. 
Recall  that  region  3  contains  the  lug  at  the  bottom  and  a  lip 
at  the  top;  so  some  metal  flow  is  bound  to  occur  to  fill 
these  undefined  features.  Regions  1  and  2  have  clearances 
from  the  die  walls  and  lateral  metal  flow  is  present  during 
forging. 


117 


Table  3«  Volumes  of  regions  in  the  weapon  component. 


Total  measured  volume  of  part  =  21.52  cc  =  1.3132  in3 

Volume  of  central  hole  =  0.0505  i n 3 

Total  volume  of  finished  part  =  1.3132  +  0.0505 

=  1.3637  in3 


Computed  volumes  for  the  described  part: 
Zone  1  Region  1 

Region  2 
Region  3 

Zone  2  Region  l 

Zone  3  Region  1 

Total  described  volume  (1/2  part) 


=  0.2202  in3 
=  0.1568  in3 
«=  0.1622  in3 
=  0.04729  in3 
=  0.06294  in3 

=  0.6494  in3 


Described  volume  of  finished  part 
Volume  of  features  not  described 


Volume  to  be  added  to  Zone  1 


=  1.2988  in3 

-  1.3637  -  1.2988 
«  0.0649  in3 

-  0.03245  in3 
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bince  each  region  of  the  preform  has  the  required  mass 
of  the  finished  shape,  densif ication  problems  are  avoided. 

Die  wear  at  the  corners  is  avoided  as  there  is  no  substantial 
flow  around  them.  Free  surface  fracture  occurrence  is  also 
inhibited  due  to  the  presence  of  compression  over  the  entire 
volume  of  the  preform.  Hence,  this  is  a  desirable  preform. 

The  clearances  for  all  the  regions  in  the  part  are 
computed  for  a  stroke  of  0.5  inches.  This  was  explicitly 
chosen  to  allow  for  a  small  clearance  in  region  3  to  permit 
the  preform  to  fall  into  the  dies  before  forging.  The 
assembled  preform  is  shown  in  Figure  35.  The  broken  lines 
indicate  the  die  walls. 

The  preform  shown  in  Figure  35  will  lead  to  lap  defects 
u j<  to  the  sharp  transition  in  the  profile  from  one  region  to 
another.  Generous  blending  is  required  to  prevent  this 
uefect.  This  is  done  externally.  The  blended  preform  is 
shown.  Figure  36,  in  isometric.  The  overall  metal  flow 
curing  forging  is  shown  in  Figure  37.  The  dimensioned 
preform  drawing  is  given  in  Figure  38.  This  is  the  preform 
that  is  trial  forged. 

4.8.3  Prototype  Forging 

Six  preforms  were  machined  according  to  the  preform 
arawing ,  Figure  38.  The  material  was  4643  steel  powder, 
compacted  to  80%  density  and  sintered  at  2050F.  Lach 
machined  preform  was  coated  with  graphite,  reheated  to  1800F 
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Figure  35-  Assembled  preform  (unblended)  showing  the  respective 
cl earances . 


FORGING  DIRECTION 
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Figure  37. 


Metal  flow  directions  in  the  designed  preform  during  forging. 


122 


i 

8*0  =i 

E 

a 

Eh 

M 

S 

s 

cn 

cx 

o 

w 

o 

Q 

Px, 

w 

§ 

s 

§ 

a, 

M 

o 

cc 

Ej 

P*-» 

M 

K 

< 

s 

PL, 

O 

Figure  3^.  Dimensioned  preform  drawing  for  prototype  Gorging  of  line  part. 
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and  forged  in  a  hydraulic  press.  Prior  to  forging , graphite 
in  water  was  sprayed  on  the  dies  for  lubrication. 

The  forging  trials  were  successful .  ( ^ ^  There  were  no 
defects  in  the  forgings.  The  lugs,  lip  etc.  were  fully 
formed.  Macrographs  taken  at  various  sections  indicated  good 
densif ication,  except  for  some  residual  porosity  at  the 
surfaces . 

A  flat  preform  was  also  forged  to  ascertain  the  defect 
formation.  As  expected,  it  resulted  in  considerable 
fractures  and  lap  defects.  This  experiment  served  to  confirm 
the  importance  of  preform  design  in  powder  forging. 


124 


5.0  CONCLUSIONS  AND  RECOMMENDATIONS 

A  computer-aided  preform  design  technique  has  been 
developed  and  applied  to  obtain  a  preform  for  a  weapon 
component.  The  designed  preform  was  forged  to  a  fully 
densified,  defect  free  part.  The  goals,  set  out  in  the 
project  have  been  achieved.  Specifically,  the  goals  were  to 
develop  a  general  approach,  as  far  as  possible,  to  input  the 
data  regarding  the  shape  of  the  finished  part  and  material; 
develop  a  preform  design  scheme  for  the  part;  design  the 
preform  for  a  specific  part  using  the  CAD  approach;  forge 
the  preform  designed;  and  critique  the  merit  and  limitations 
of  this  approach. 

The  CAD  program  developed  consists  of  two  phases.  The 
first  phase  is  the  description  of  the  part  geometry,  by  a 
simple  but  efficient  method.  Currently  available  techniques 
for  geometry  description  include  sophisticated  drafting 
packages  (IBM/Lockheed  CADAM)  and  the  NC  part  description 
program  APT.  The  use  of  these  were  rejected  due  to  the  con¬ 
siderable  training  required  for  the  potential  user.  A  part 
description  scheme,  called  PADEL,  was  developed  for  this 
project.  PADEL  has  features  which  are  suited  to  forging 
applications.  Specifically,  forging  direction  specification, 
sub-division  to  regions  and  the  identification  of  punch/die 
profiles.  It  is  relatively  simple  to  use.  However,  it  has 
limitations  on  the  complexity  of  part  geometry  that  can  be 
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described.  These  limitations  are  justifiable,  as  the  goal  of 
this  research  was  not  toward  the  development  of  another 
"complete''  geometry  description  scheme. 

The  second  phase  is  the  preform  design  phase.  Preform 
design  is  accomplished  on  an  interactive  basis  where  the  user 
supplies  trial  preform  shapes  which  are  evaluated  by  the 
program  for  densif ication  and  optimum  metal  flow  modes.  By 
interactively  manipulating  the  preform  shapes,  a  desirable 
preform  is  obtained.  The  computer  evaluation  is  based  on 
empirical  guidelines  developed  through  experience  with  powder 
preform  forging.  The  interactive  approach  eliminates  most  of 
the  trial-and-error  involved  in  preform  design  as  is 
commercially  practiced. 

The  effects  of  all  process  parameters  on  metal  flow  and 
fracture  in  powder  preform  forging  is  too  complex  to  be 
incorporated  into  this  program.  As  a  result,  an  exact 
definition  of  the  preform  is  impossible.  Some  minor 
adjustments  may  be  necessary  after  a  few  forging  trials. 
Evolution  of  the  program  through  user  experience  and 
increased  understanding  of  the  mechanics  of  metal  flow  is 
expected  to  lead  to  constant  improvements. 

An  immediate  extension  that  can  be  considered  is  the 
extension  of  the  current  approach  to  axi  -  symmetric  parts.' 

The  present  program  can  also  be  used,  with  some 
modifications,  to  determine  the  punch  strokes  required  for 
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hot  repressing  of  preforms.  In  this  case,  user  shape  input 
can  be  bypassed  as  the  preform  shape  conforms  closely  to  the 
finished  shape.  A  further  extension  is  the  design  of  die 
compaction  processes,  where  extensive  shape  manipulation 
capabilities  are  not  required. 

The  present  program  is  implemented  on  a  DEC-10  computer. 
Graphics  capability  is  required,  and  the  TEKTRONIX  PLOT-10 
software  is  used  for  this  purpose.  All  programming  is  in 
FORTRAN- 10 ,  except  for  some  specialized  routines  which  are 
coded  in  assembly  language.  Some  modifications  will  be 
required  to  execute  on  a  different  computer. 

The  state  of  analytical  techniques  was  reviewed  for 
application  to  the  computer-aided  design  program.  Two 
techniques  were  developed,  the  slab  analysis  and  the  energy 
method,  in  its  simplest  form.  This  development  led  to  some 
interesting  observations.  Specifically,  the  influence  of 
friction  at  the  die/part  interface  on  the  densif ication 
behavior.  Perhaps  a  fundamental  experimental  and  analytical 
study  is  warrented  at  this  stage.  The  energy  method 
developed,  can  be  extended  using  two-dimensional  elements  to 
obtain  detailed  density  distributions  in  a  deforming  preform. 
Furthermore,  it  appears  to  be  a  promising  technique  for  the 
application  of  workability  analysis  to  preform  design. 
Considerable  research  activity  is  foreseen  in  this  area. 
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To  start  an  analytical  evaluation  of  preforms,  the 
preform  must  be  defined.  The  CAD  approach,  using  empirical 
guidelines,  can  be  used  to  specify  the  approximate  preform 
that  is  subjected  to  analytical  evaluation.  Thus  future 
activity  should  be  directed  along  two  fronts.  First,  the 
improvement  of  the  empirical  approach  to  yield  a  fairly 
reliable  preform.  Second,  the  development  of  the  energy 
method  of  a  similar  technique  capable  of  detailed  results 
towards  the  evaluation  of  fracture  and  dens  if ication  of 
preforms.  Once  these  objectives  are  accomplished,  an 
integrated  powder  forging  process  design  may  be  assembled. 
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APPENDIX  A 


Radial  Stress  Distribution  for  Ax i symmetric  Slab 

Consider  the  n-th  slab  or  zone,  as  shown  in  Figure  10. 

The  slab  is  located  between  the  radial  position  r  ,  r  ,  . 

r  n  n+1 

The  pressures  acting  on  an  infinitesmal  element  are  also 
indicated.  Taking  force  equilibrium  along  the  coordinate 
directions,  we  have 
Radial : 


(Pr+dp^_)hrd0  -  p^_(r+dr) 

-  T^r+dr/^Cosci^de  dS^^ 

-  p^ (r+dr/2)Sina^d0 
=  0 


(h+dh)d0  -  t  (r+dr/2)Cosa  d0  dS 
u  u  u 

-  p  (r+dr/2)Sina  d0  S 

u  u  u 

-  2pQ (h+dh/2) Sin (d0/2)  dr 


Assuming  that  Sin  d0/2  =  d0/2  and  Pg  =  Pr  /  neglecting 
higher  order  differentials  the  above  can  be  simplified  to 


dPr/dr  =  {pr  dh/dr  +  +  (p^ano^+p^tana  )  }  /  h 


(A-l) 
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Axial : 

(i)  Upper  Interface. 

P  (r+dr/2)d0  dr  =  p  (r+dr/2)Cosa  dO  dS  -  x  (r+dr/2) Sina  d0  dS 
z  u  u  u  u  u  u 

which  simplifies  to 


P 


u 


p  +  x  tana 
z  u  u 


(ii)  Lower  interface. 


P1  = 


+  T1tano1 


Now  by  the  geometry  of  the  element, 
-dh/dr  =  (tana^  +  tana^) 


El  iminating  PU»P1  and  dh/dr  from  (A-l),  and  changing 

pressures  to  stresses  according  to  a  =  -p  ,  a  =  -p  we  have 

r  r  z  z 

da  /dr  =  -  { (a  -a  ) (tana  +tana  )  +  (x  Sec2a  +  x  Sec2a  ) }  /  h 
r  ‘  iZ  UJ.  null 


with 


r  =  r  ;  (a  )  =  (a  ) 

n+1  r  n+1  r  r=r 


(A-2) 


n+1 


The  interface  shear  is  represented  by  a  friction  factor 
according  to 

x  =  mk  ( A  —  3  ) 

where  k  :Shear  yield  of  the  porous  compact, 
m  :  friction  factor. 
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APPENDIX  B 


Rectangular  Element  Equations  for  Energy  Method 


Consider  the  m-th  element  flanked  by  sides  i  and  j 


Then , 

(i)  Velocity  vector  interpolation. 


T 

I*  ,  =  G  V.  . 
(m)  ~  (m) 


(B-l) 


where 


T 

Y(m) 


v. 
-  1 


V 


GT  =  Lq, 


=  (l+s)/2  j  '  q  *  (l-s)/2  ;  -1  <  s  <  1 


(ii)  Coordinate  Transformation. 


r  =  q.r.  +  q.r. 
i  i  3  3 


(B-2) 


(iii)  Strain-rate  vector. 


'  (m)  ~  (m) 


3u/3r 

u/r 


(B-3) 


where 


B  = 


_1/rji  1/rji 


Vr  yr 
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(iv)  Equivalent  Strain-rate . 
Simplifying  (3-25)  we  get 


e2 

~eq (m) 


•  ip  • 

2  e  Ae  + 


•  • 

e  Ce  +  2v 


T* 

c  e 


3 (1+v) 
1-v 


3 (l-2v) 


(1+v) (l-2v) 


(1+v) (l-2v) 


£2 

z 


where 

A 


1.0  -0.5 

-0.5  1.0 


c  = 


1.0 


1.0 


1.0  1.0 


c  = 


1.0  1.0 


Substituting  (B-3)  and  simplifying 


E2  =  v  Xv  +  2a,  c  Bv  +  a  (B-4) 

eq  ~  1  0 

where  the  element  subscript  (m)  is  dropped  and 


x  =  2  btab  +  l  btcb 

3 (1+v)  3 (l-2v) 


a  =  _ v _  e 

(1+v) (l-2v)  2 


•  o 

a  =  1-v _  £ 

(1+v) (l-2v)  2 


(v)  Element  level  functional . 

The  functional  at  the  element  level  is 


/  o  e  dv 
v  eq  eq 


f 


T 

I  G  v  dS 


S 


F 
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Substituting  (B-4),  we  have 

=  f  oeq(  yTXv  +  2a1cTBy  +  aQ)*  dV  -  I  xGTy  dS 
V  SF 

and  hence, 

-/  o  (  Xy  +  a  BTc)/(e2  )*  dv  -  /  tG  dS 
v  «I  1  eq 


(B-5) 


(vi)  Linearization . 

At  the  n-th  iteration  on  velocity  components,  suppose 


' (n)  ■  *<„-!>  +  4V(„) 


(B-6) 


Substituting  (B-6)  into  (B-4)  and  retaining  first  order  terms 

ln  AY(n)' 


e2 

eq(n) 


•  2 

£  ^ 

eq (n-1) 


T  T  T 

+  2(  v  X  +  Q  )  Ay 


(n) 


e2 


eq(n-l)  +  2a(n-l)  ^ (n) 


where 


Q  =  a  B  c 

1 


and 


a .  .  =  Xv  ,.  +  o 

(n-1)  -  (n-1)  * 


Hence , 

•-1 

e 

eq 


u2  r4 

eq 


.-1  * -3  T 

'eq  £eq  3 (n-1)  ~(n) 


(B-7) 
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The  volume  integrand  of  ( B— 5 )  at  the  n-th  iteration  is 


=  a  e_1 ,  .  (  Xv,  .  +  Q  ) 

eq  eq(n)  (n) 


Substituting  (13-6)  and  ( B— 7 )  we  have, 


°eq  Eeq(n-1)  a(n-l)  +  Ceq(n-1)  a(n-l)a(n-l)  ^~(n)  ^ 


Defining  r  =/  tG  dS ,  (B-5)  becomes 
S„ 


9<J>/9V  =  a  P.  Av.  _.  +  a  H.  ..  -  T 

-  eq  (n-1)  ~ (n-1)  eq  (n-1) 


(B-8) 


where 


(n-1) 


f  e 


-1 

eq (n-1) 


{  X  -  e 


-2  T 

eq(n-l)  3 (n-1)  a(n-l) 


}  dV 


H/  i \  =  /  e  1 ,  , .  a .  . .  dV 

(n-1)  eq(n-l)  (n-1) 


Assembling  (B-8)  we  have 

M 

E  d<f>/  9v  =0 

m=l 

The  integrals  are  numerically  evaluated  using  a  2-point 
Gaussian  quadrature.  The  coefficient  matrix  of  the  system  of 
linear  eouations  is  tri-diagonal  and  is  solved  using 
recursion  relationships. 
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PADEL  STATEMENT  FORMATS 

C.l  Part  description  space  statements 

(i)  Point  statements . 

P POINT  /  X,  Y,  Z 

PPOINT  /  plane  1,  plane  2,  plane  _3 


(ii)  Plane  statements . 

PLANE  /  ppoint ,  PARLEL ,  plane 

PLANE  /  ppoint  1,  ppoint  2,  PERPTG,  XYPLN 

PLANE  /  PARLEL,  plane,  [ X LARGE ,  XSMALL ,  YLARGE , 

Y LARGE ,  ZLARGE,  Z3MALL]  ,  d 

The  keyword  "CUTPLN"  may  be  substituted  for  "PLANE"  to 
define  cut  planes.  Note  that  the  positive  unit  normals  for 
the  relative  plane  definitions  are  in  the  same  sense  as  the 
parent  plane.  For  the  second  definition,  the  sense  of  the 
normal  is  defined  by 

n  x  v  =  k 

where,  n  =  positive  unit  normal. 

v  =  unit  vector  from  ppoint  1  to  ppoint  2  projected 
on  XYPLN. 

k  =  positive  unit  normal  of  XYPLN. 


C.2  Primary  definition  plane  statements 


(i)  Point  statements. 


POINT 

/ 

>1 

X 

POINT 

/ 

INTOF,  line  1,  line  2 

POINT 

/ 

( XLARGE ,  XSMALL , 

YLARGE, 

YSMALL],  INTOF, 

line,  circle 

POINT 

/ 

[XLARGE,  XSMALL, 

YLARGE, 

YSMALL],  INTOF, 

circle  1,  circle 

2 

(ii)  Line  statements . 


LINE 

/ 

x,  y,  x,  y 

LINE 

/ 

point  1,  point  2 

LINE 

/ 

[XAXIS,  YAXIS) 

LINE 

/ 

point ,  [LEFT,  RIGHT],  TANTO, 

LINE 

/ 

point,  [PARLEL,  PERPTO] ,  line 

LINE 

/ 

point,  ATANGL,  ©  ,line 

LINE 

/ 

PARLEL,  line,  [XLARGE,  XSMALL 

YLARGE,  YSMALL] ,  d 
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A 


( i i i )  Circle  statements . 

CIRCLE  /  x,  y,  r 

CIRCLE  /  point  1,  point  2,  point  3 
CIRCLE  /  CENTER,  Point,  RADIUS,  r 
CIRCLE  /  CENTER,  X,  y,  RADIUS,  r 
CIRCLE  /  CENTER,  point  1,  point  2 
CIRCLE  /  CENTER,  point,  TANTO,  line 
CIRCLE  /  CENTER,  point,  [LARGE,  SMALL] 
TANTO,  circle 

CIRCLE  /  TANTO,  line ,  [XLARGE ,  XSMALL , 
YSMALL] ,  point,  RADIUS,  r 


C.3  Non-definition  statements 

PARTNAME/  (user  defined  name> 
FINISH/ 

ZONE/n 

ENDZON/n 

SYMMETRY/  plane  1.,  plane  2,  plane  3 

BETWEEN/  plane  1,  plane  2 

FRAME/  plane  _1,  plane  2 

REVERSE/  plane  1.,  plane  2,  plane  3 

INCLUDE/  plane  1.,  plane  2 

TOP/  point 

BCT/  point 

REGION/  line  1,  line  2,... 


YLARGE , 
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GOFKD/  line  1,  line  2,... 

GOFWD/  line  1,  line  2  , . . .circle 
GOCLK/  circle  1,  circle  2 
GOCLK/  circle  _1,  line 
GOACLK/  circle  1,  circle  2 
GOACLK/  circle  1^  line 
ADD/ 
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APPENDIX  D 


PADEL  SOURCE  STATEMENTS  FOR  A  PART 


;*****  TYPICAL  SOURCE  FILE  FOR  PADEL  ***** 

m 

9 

; PART  TITLE  (MAX  30  CHARACTERS) 

f 

PARTNAME/WEAPON-COMPONENT 

9 

; ORIENT  PART  IN  DESCRIPTION  SPACE. 

; SPECIFY  THE  DEGREE  OF  SYMMETRY;  IF  ANY. 

• 

9 

SYMMETRY/YZPLN 

t 

; SUP-DIVISION  INTO  ZONES. 

; DEFINE  APPROPRIATE  CUT  PLANES. 

9 

CTPLl =  CUTPLN/PARLEL,ZXPLN,YLARGE,0.G16 
CTPL2=  CUTPLN/PARLEL , ZXPLN,Y LARGE, 0.241 

THE  CUT  PLANES  ARE  DEFINED  WITH  RESPECT  TO  THE 
PART  COORDINATE  SYSTEM.  NOW  FOP  EACH  INTENDED  ZONE, 
IDENTIFY  THE  FOLLOWING: 

(1)  PRIMARY  AND  SECONDARY  DEFINITION  PLANES. 

(2)  LOCAL  COORDINATE  PLANES. 

THE  COORDINATE  PLANES  UPON  INTERSECTION  WITH  THE 
PRIMARY  DEFINITION  PLANES  DEFINE  THE  X-AXIS  AND 
Y-AXIS  FOR  THE  ZONE  CONTOURS.  IT  IS  IMPERATIVE 
THAT  THE  PROGRAMMER  KNOW  THE  SENSE  OF  THE  POSITIVE 
UNIT  NORMALS  OF  THE  PLANES  SINCE  THEY  REPRESENT 
THE  POSITIVE  AXES  OF  THE  ZONE  FRAME. 

THE  PERTINENT  PLANES  FOR  THE  ZONES  ARE: 

ZONE  1. 

PRIMARY  PLANE. . .PRIM1 

P  R I M 1 =  PLANE/PARLEL,YZPLN,XLARGE,0 .515 

SECONDARY  PLANE  ... YZ PLN 
LOCAL  X-AXIS  PLANE. . .CTPLl 
LOCAL  Y-AXIS  PLANE... LOCY 

LOC Y=  PLAN E/PA RLEL,XYPLN,X LARGE, 4 .0 
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ZONE  2  . 

PRIMARY  PLANE. . .PRIM2 

PPIM2=  PLANE/ PARE EL, YZPLN ,X LARGE , 0 .215 

SECONDARY  PLANE ... YZPLN 
LOCAL  X-AXIS  PLANE. . .CTPL2 
LOCAL  Y-AXIS  PLANE... LOCY 


ZONE  3. 

PRIMARY  PLANE. .. ZXPLN  (REVERSED) 

SECONDARY  PLANE ... CTPL2 
LOCAL  X-AXIS  PLANE. . .YZPLN 
LOCAL  Y-AXIS  PLANE  ... LOCY3 

LOC Y3=  PLANE/PARLEL#  LOCY , ZSMALL ,0.563 


END  OF  ALL  PLANE  DEFINITIONS. 


START  DEFINITION  OF  ZONES. 

ZONE/I 

INDICATE  THE  DEFINITION  PLANES. 

BETWEEN/PRIMl , YZPLN 
INDICATE  THE  COORDINATE  PLANES. 

FRAME/CTPLl , LOCY 
INDICATE  INTERFACE  PLANES. 

INCLUDE/ CTPLl 

START  OF  CONTOUR  DEFINITION. 

PT0  =  POINT/INTOF,XAXIS,YAXIS 
LNl=  LINE/PARLE L  ,  YAXIS ,XLARGE , 0 .663 
PTCI=  POINT/0.798,-0.375 
PTC2=  POINT/1.798,-0.375 
CIRl=  CIRCLE/CENTER, PTC1 , RADIUS, 0.1 35 
CIR2=  CIRCLE/CENTER, PTC2 , RADIUS ,0 . 135 
LN  2=  LINE/PARLEL,XAXIS,YSMALL,0.51 
LN3=  LINE/PARLEL , LN 1 , XLARGE , 1 . 25 
PTl=  POINT/2.268,-0.935 


LN4  =  LIN-E/PTl  ,  PARLEL , YAXI S 
LN5=  LINE/PT1 , ATANGL ,-l 5 , XAX IS 
PT2=  POINT/0.0,-0.61 
LN6  =  LINE/PT2 , ATANGL, -6 ,XAXIS 

CONTOUR  SEQUENCE  SPECIFICATION. 

TOP/PT0 

GOFWD/XAX IS  ,  LN1  ,  C I  Pi 
GOACLK/CIRl ,  LN2 
GOFWD/LN2 ,  CIR2 
GOACLK/CIR2 ,LN3 
GOFWD/LN3  , XAXI5 , LN4 

BOT/PT2 

GOFWD/LN6 , LN5 , LN4 
REGION  DEFINITION  FOR  ZONE 
REGION/LNl ,  LN3 

DOES  ANY  MASS  HAVE  TO  EE  USER  INPUT?  IF  YES  THEN 
ADD/ 

ALL  INPUT  FOP  ZONE  I  COMPLETE  :  TERMINATE. 
ENDZON/1 


INPUT  FOR  ZONE  II. 

ZONE/2 

BETWEEN/PRIM2 , YZPLN 
FRAME/CTPL2 , LOCY 
INCLUDE/CTPL1 ,CTPL2 

PT0*  POINT/0.0,0.0 
PTl=  POINT/0.375,0.0 
PT3=  POINT/0.375,-0.61 
PT4=  POINT/0.0,-0.563 

TOP/PT0 

GOTO/PTl 

BOT/PT4 

GOTO/PT3 


ENDZON/2 


IN' PU1  IOP  ZONK  III 


ihi 


ZONE/ 3 

BETViEEN/ZXPLN  , CTPL2 
REVERSE/ZXPLN 
FRAML/YZPLN , LOCY3 
INCLUDE/YZPLN 

PT0®  POINT/0.0,0.0 

PT1=  POINT/ 0.499, 0.563 

LNl=  LINE/PH  ,PARLEL,  YAXIS 

PT2*  POINT/0.16,0.0 

LN  2*  LINE/PT2 , ATANGL, 19.0,XAXIS 

PT3*  POI NT/I NTOF  , LN2 , LN 1 

PT4  =  POINT/0.0,0.563 


TOP/FT4 
GOTO/ PTl 

DOT/PT0 

GOTO/PT2 

COTO/PT3 

ENDZON/3 


D I ,S C  R 1 PT I  ON  COMPLETE  . 


FINISH/ 


1JU2 
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